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Scilab numbering policy used in this document and the relation to the 
above book. 

Exa Example (Solved example) 

Eqn Equation (Particular equation of the above book) 

AP Appendix to Example(Scilab Code that is an Appednix to a particular 
Example of the above book) 

For example, Exa 3.51 means solved example 3.51 of this book. Sec 2.3 means 
a sc-ilab code whose theory is explained in Section 2.3 of the book. 
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Chapter 5 

Ideal Cycles and their Analysis 


Scilab code Exa 5.1 Calculation of MEP and Efficiency 
clc ; 

p 1 = 1 ; // Pressure before compression in bar 

Tl=350; // Temperature before compression in kelvin 

T3 = 2000; // Temperature after combustion in kelvin 

rp=1.3; // Pressure ratio 

Cp = 1.005; // Specific heat at constant pressure in 
kJ/kg K 

r = 1.4; // Specific heat ratio 

R = 287; // Characteristic gas constant in J/kg K 

T2 = Tl*(rp) “ ( (r-1)/r) ; // Temperature at the end of 
t h e c o m p r e s s i o n 

T4 = T3* (1/rp) ~ ((r-1)/r) ; // Temperature after 
e x p a n s i o n 

Wc = Cp* (T2-T1 ) ; // Work done during compression 

WT = Cp* (T3-T4) ; // Work done during expansion 

WN = WT-Wc; // Net work done 

p2 = rp*pl ; // Pressure at state 2 

p3 = p2; p4 = pl ; // Constant pressure process 

VI =R*T1 / (pl * 10~5) ; // specific Volume at state 1 

V2 = R*T2 / (p2 * 10 ~ 5) ; // specific Yolume at state 2 
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V3 = R*T3 / ( p3 * 10 ~ 5) ; // specific Volume at state 3 
V4 = R*T4/(p4*10~5) ; // specific Volume at state 4 
imep = WN* 10 ~ 3/(V4 - V2 ); // Mean effective pressure 
q=Cp * ( T3-T2 ) ; // Heat supplied 
eff=WN/q; // Efficiency of a Jonle cycle 
disp (” bar”,imep* 10“-5 , ” Mean effective pressnre = ”) 
> 

disp (”%” , ef f *100 , ” Efficiency of a Joule cycle =”); 


Scilab code Exa 5.2 Calculation of Improvement in Efficiency 
clc ; 

p 1 = 1 ; // Pressure before compression in bar 

Tl=350; // Temperature before compression in kelvin 

T3 = 2000; // Temperature after combustion in kelvin 

rp=1.3; // Pressure ratio 

Cp = 1.005; // Specific lieat at constant pressure in 
kJ/kg K 

r = 1.4; // Specific heat ratio 

R = 287; // Characteristic gas constant in J/kg K 

T2 = T1 * (rp) “ ((r-1)/r) ; // Temperature at the end of 
t h e c o m p r e s s i o n 

T4 = T3* (1/rp) “ ( (r-1)/r) ; // Temperature after 
e x p a n s i o n 

Wc = Cp* (T2-T1 ) ; // Work done during compression 
WT = Cp* (T3-T4) ; // Work done during expansion 
WN = WT-Wc; // Net work done 
T5=T4; // For a perfect heat exchange 

q=Cp * (T3-T5 ); // Heat added 

eff2 = WN/q; // Efficiency of a modified Joule cycle 
ef f 1 = 0.072220534; // Efficiency of a joule cycle 

disp (”%”,eff2*100,”Efficiency of a modified Joule 
cyclc = ” ) ; 

disp ( ef f 2/ef f 1 , ” Improvement in efficiency = ”); 
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Scilab code Exa 5.3 calculation of net power output of the cycle 
clc ; 

rp = 6; // Pressure ratio 

Tl=300; // Inlet air temperature to the compressor 

i n k e 1 v i n 

T3 = 577 + 273; // Inlet temperature of air at turbine 

i n k e 1 v i n 

Vr = 240; // Volume rate in m'3/s 

Cp = 1.005; // Specific heat at constant pressure in 
kJ/kg K 

r = 1.4; // Specific heat ratio 

R = 287; // Characteristic gas constant in J/kg K 

p1=1; // pressure at state 1 in bar 

T2 = Tl*(rp) ~ ((r-1)/r) ; // Temperature at the end of 
t h e c o m p r e s s i o n 

T4 = T3* (1/rp) “ ((r-1)/r) ; // Temperature after 
expansion 

Wc = Cp* (T2-T1 ) ; // Work done during compression 
WT = Cp* (T3-T4) ; // Work done during expansion 
WN = WT-Wc; // Net work done 
q=Cp * (T3-T2 ) ; // Heat supplied 

rowl=pl*10~5/(R*Tl); // Density of air at state 1 
P=WN*Vr*rowl ; // Power output 
eff=WN/q; // Efficiency of a cycle 

disp (”MW (roundoff error ) ” , P/1000 , ” Power Output = 
”); 

disp (”%” , ef f *100 , ” Efficiency of a cyclc = ”); 


Scilab code Exa 5.4 Calculation of Ehiciency and work of compression 
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1 clc ; 

2 Tl=300; // Inlet air temperature to the compressor 

i n k e 1 v i n 

3 pl=l; // pressure at state 1 in bar 

4 T2=475; // Temperature at discharge in kelvin 

5 p2 =5; // Pressure at state 2 

6 T5 = 655; // Temperature after heat exchanger in 

k e 1 v i n 

7 T3 = 870 + 273; // Temperature at he turbine inlet in 

k e 1 v i n 

8 T4 = 450 + 273; // Temperature after turbine in kelvin 

9 Cp = 1.005; // Specific lieat at constant pressure in 

kJ/kg K 

10 r = 1.4; // Specific lieat ratio 

11 R = 287; // Character ist ic gas constant in J/kg K 

12 

13 Wc = Cp* (T2-T1) ; // Work done during compression 

14 WT = Cp* (T3-T4) ; // Work done during expansion 

15 WN = WT-Wc; // Net work done 

16 q=Cp* (T3-T5) ; // Heat supplied 

17 eff=WN/q; // Efficiency of a cycle 

18 

19 disp (”kJ/kg” ,WN ,” ( i ) . The output per kg of air = ”) 

) 

20 disp (”%” ,eff *100,” ( ii ) . The efficiency of the cycle 

= ”) ; 

21 disp (” kJ/kg” , Wc , ” ( i i i ) . The work required to drive 

the compressor = ”); 


Scilab code Exa 5.5 Calculation of Thermal Efficiency 

1 clc ; 

2 p1=1.4; // Pressure at state 1 in bar 

3 Tl=310; // Temperature at state 1 in kelvin 

4 rp=5; // Pressure ratio 
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Traax = 1050; // Maximum temperatuer in kelvin 

WN = 3000; // Net output in kW 

Cp = 1.005; // Specific lieat at constant pressure in 
kJ/kg K 

r = 1.4; // Specific lieat ratio 

R = 287; // Characteristic gas constant in J/kg K 

T3=Tmax; 

T2 = Tl*(rp) “ ((r-1)/r) ; // Temperature at the state 2 
T4 = T3/ (rp) “ ((r-1)/r) ; // Temperature at the state 4 
T5=T4; // As regenerator effectiveness in 100 % 
m = WN/( Cp * ( (T3-T4) - (T2 -T1) ) ) ; // mass flow rate of 
air 

ef f = (T3-T4-T2 + T1) / (T3-T5) ; // Efficiency of a cycle 
disp (”%” , ef f *100 , ” ( i ) . Thermal efficiency of the 

cycle = ” ) ; 

disp (”kg/min (roundoff error)”,m*60,”(ii). The 
mass flow rate of air per minute = ”); 


Scilab code Exa 5.6 Calculation of Pressure ratio of compressor and tur- 
bine 

clc ; 

Tl=290; // Compressor inlet temperature in kelvin 

T2=460; // Compressor outlet temperature in kelvin 

T3 = 900 + 273; // Turbine inlet temperature in kelvin 

T4 = 467 + 273; // Turbine outlet temperature in kelvin 

Cp = 1.005; // Specific heat at constant pressure in 

kJ/kg K 

r = 1.4; // Specific heat ratio 

R = 287; // Characteristic gas constant in J/kg K 

c = T2/Tl ; // Temperature ratio 

rpc = c ~ ( r / ( r - 1 ) ) ; // Compression ratio 

WN = (Cp*((T3-T4)-(T2-T1) ) ) ; // Specific power 
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13 T5 = T4; // Assuming regenerator e f f ec t i v en e s s to be 

100 % 

14 ef f =WN/ ( Cp * ( T3-T5 ) ) ; // Overall efficiency of the 

cycle 

15 Wc = Cp*(T2-Tl); // Work required to drive the 

compressor 

16 rpt = (T3/T4) “(r/(r-l) ) ; // Turbine pressure ratio 

17 disp (rpt,” Turbine pressure ratio = ”,rpc,” 

Compressor pressure ratio = ”,”(i).”); 

18 disp (” kJ/kg” , WN , ” ( i i ) . Specific power output = ”); 

19 disp (”%” ,eff *100 , ”(iii). Overall efficiency of the 

cycle = ” ) ; 

20 disp (” kJ/kg” , Wc , ” (iv). Work required to drive the 

compressor = ” ) ; 


Scilab code Exa 5.7 Calculation of temperature drop across the turbine 

1 clc ; 

2 nW_WT = 0.563; // Ratio of net work to turbine work 

3 Tl=300; // Inlet temperature to the compressor in 

k e 1 v i n 

4 ef f =0 . 35 ; // Thermal efficiency of the unit 

5 m = 10; // massflow rate in kg/s 

6 Cp = l; // Specific heat at constant pressure in kJ/kg 

K 

7 r = 1.4; // Specific lieat ratio 

8 

9 c = l/(l-eff); // For ideal simple cycle 

10 T2 = Tl*c; // Temperature at state 2 

11 Wc = Cp* (T2-T1) ; // Compressor work 

12 WT = Wc / ( 1 - nW_WT ) ; // Turbine work 

13 WN = WT-Wc; // Net work 

14 q=WN/eff; // Net lieat supplied per kg of air 

15 T3 = (q/Cp) +T2 ; // Temperature at state 3 

16 T4 = T3/c; // Temperature at state 4 
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17 T3_T4 = T3-T4; // Temperature drop across the turbine 

18 disp (”K” , T3_T4 , ” Temperature drop across the turbine 

= ”); 


Scilab code Exa 5.8 Calculation of turbine pressure ratio 

1 clc ; 

2 p=336.5; //specific power output of a turbine in kW/ 

kg 

3 T4 = 700; // Temperature at turbine outlet in kelvin 

4 Cp = l; // Specific heat at constant pressure in kJ/kg 

K 

5 Cv = 0.717; // Specific lieat at constant volume in kJ/ 

kg K 

6 

7 r = Cp/Cv; // Specific heat ratio 

8 T3 = T4 + (p/Cp) ; // Temperature at turbine inlet 

9 p3_p4=(T3/T4) ~ (r/(r-l)); // Pressure ratio across 

the turbine 

10 disp (round (p3_p4), ”Pressure ratio across the 
turbine = ”); 


Scilab code Exa 5.9 Estimation of thermal efficiency of the plant 

1 clc ; 

2 Tl=300; // Minimum operating temperature in kelvin 

3 T3 = 900; // Maximum operating temperature in kelvin 

4 p 1 = 1; // Minimum pressure in bar 

5 p3=4; // Maximum pressure in bar 

6 m = 1600; // Mass flowrate in kg/min 

7 r = 1.4; // Specific heat ratio 

8 Cp = 1.005; // Specific lieat at constant pressure in 

kJ/kg K 
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p2=p3; p4=pl; // Constant pressure process 

c = (p2/pl)~((r-l)/r) ; 

eff = (l-l/c); // The efficiency of the cycle 
t = T3/Tl; // ratio of maximum and minimum temperature 
W = Cp*Tl * (t * (1-1/c) - (c-1) ) ; // Work output per kg of 
air 

P=(m/60)*W; // Shaft power available 
disp (”%” , ef f *100 , ” Thermal efficiency of the plant 
= ”); 

disp (”kW (roundoff err or ) ” , P , ” Shaft power 
available for external Load = ”); 


Scilab code Exa 5.10 Determination of the cycle thermal efhciency 
clc ; 

Tl = 15 + 273; // Inlet temperature of air at compressor 

i n 1 e t i n k e 1 v i n 

rp = 6; // Compressor pressure ratio 

T3=750 + 273; // Maximum permissible temperature in 

k e 1 v i n 

T5 = T3; // After reheat 

Cp = 1.005; // Specific lieat at constant pressure in 
kJ/kg K 

r = 1.4; // Specific lieat ratio 

c = rp~((r-1)/r) ; 

T2 = Tl*c; // Temperature at state 2 

p3_p4 = sqrt (rp); // For maximum expansion work 

T4 = T3/(p3_p4) ~ ((r-1)/r) ; // Temperature at state 4 

T6 = T4; // As pressure ratio is same 

Wc = Cp* (T2-T1 ) ; // Compressor work 

WT = Cp*(T3-T4)+Cp*(T5-T6) ; // Turbine work 

T7=T4; // Because of 100% regeneration 

q=Cp * ( T3-T7) +Cp * ( T5-T4 ) ; // Heat supplied 
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18 WN = WT-Wc; // Net work done 

19 eff=WN/q; // Efficiency of the plant 

20 Wratio=WN/WT; // Work ratio 

21 disp (”kJ/kg of air”,q,”Heat supplied = ”); 

22 disp (”kW (roundoff er r or ) ” , WN , ” Net shaft work = ” 

) ; 

23 disp (”%” ,eff*100,”The cycle thermal efficiency = ”) 

> 

24 disp (Wratio , ” Work ratio = ”); 


Scilab code Exa 5.11 Calculation of Efhciency under conditions giving 
maximum work 

1 clc ; 

2 Tmin = 5 + 273; // Minimum operating temperature in 

k e 1 v i n 

3 Tmax = 839 + 273; // Maximum operating temperature in 

k e 1 v i n 

4 Cp = 1.005; // Specific heat at constant pressure in 

kJ/kg K 

5 r = 1.4; // Specific heat ratio 

6 

7 ef f _carnot = 1 -Tmin/Tmax ; // Efficiency of the carnot 

cycle 

8 c=1/(1-eff_carnot); 

9 p2_pl=c“(r/(r-l)); // Pressure ratio 

10 disp (p2_pl , ” ( i ) . Pressure ratio at wliich efficiency 

equals Carnot cycle efficiency = ”); 

11 t = Tmax/Tmin; // Temperature ratio 

12 // Pressure ratio for maximum work is obtained when 

13 c = sqrt (t ) ; 

14 p2_pl = c~(r/(r-l)) ; // Pressure ratio 

15 eff = l-l/c;// Efficiency at maximum work output 

16 disp (p2_pl , ” ( i i ) . Pressure ratio at wliicli maximum 

work is obtained = ”); 
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disp (”%” , ef f *100 , ” ( i i i ) . Efficiency at maximum work 
output = ” ) ; 


Scilab code Exa 5.12 Comparison of basic cycle with moclihed cycles 
clc ; 

rp=4;// Overall pressure ratio 
Tl=300; // Temperature at state 1 in kelvin 

T3 = 1000; // Temperature at state 3 in kelvin 

Cp = l; // Specific heat at constant pressure in kJ/kg 
K 

Cv = 0.717; // Specific lieat at constant volume in kJ/ 

kg K 

/ / Basic cycle 

r = Cp/Cv; // Specific lieat ratio 
c=rp~((r-1)/r); 
t = T3/Tl; // Temperature ratio 

WN = Cp*Tl * (t * (1-1/ c) - (c-1) ) ; // Net work output 
ef f = (1-1/c) * 100 ; // Efficiency of the cycle 

// Basic cycle with heat exchanger 

WN_he = WN; 

ef f _he = ( 1 - c/t ) * 100 ; // Efficiency of the cycle with 

heat exchanger 

dev_WN1 = ( WN_he-WN )* 100/WN ; //Percentage deviation of 
Net work from basic cycle 
dev_ef f 1 = ( ef f _he - ef f ) * 100/ef f ; // Percentage 
deviation of efficiency from basic cycle 

// Basic cycle with intercooled compressor 

WN_ic = (Cp*Tl)*(t*(l-l/c) -2*(sqrt (c)-1)) ; 
eff_ic=(l-(((t/c)+sqrt(c)-2)/(t-sqrt(c))))*100; 
dev_WN2 = (WN_ic-WN)* 100/WN; //Percentage deviation of 
Net work from basic cycle 
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25 dev_eff 2 = ( eff_ic - eff )* 100/eff ; // Percentage 

devi at ion of efficiency from basic cycle 

26 

27 // Basic cycle with heat exchanger and intercooled 

compressor 

28 WN_iche=WN_ic ; 

29 eff_iche=(l-((2*( sqrt (c)-l))/(t*(l-l/c))))*100; 

30 dev_WN3 = (WN_ i che - WN) * 100/WN ; //Percentage deviation 

of Net work from basic cycle 

31 dev_eff3 = (eff_iche-eff)* 100/eff; // Percentage 

deviation of efficiency from basic cycle 

32 

33 printf (”Cycle \t\t\t\t\t\t WN(kJ/kg) \t\ 

tefficiency (in percentage)\t\t percentage 
Change in WN \ t \ t per cent age change in efficiency” 
); 

34 printf (” \n\t\t\t\t\t\t ( roundoff error) \t( 

roundoff error) \t\t\t (roundoff error)\t\t\t\t 
(roundoff error )” ) ; 

35 printf (”\n\nBasci cycle \t\t\t\t\t %f \t\t\t %f\t\t 

\t\t\t — \t\t\t\t\t — ”,WN,eff); 

36 printf (”\n\nWith Heat Exchanger \t\t\t\t %f \t\t\t 

%f\t\t\t\t\t %f \t\t\t %f” , WN_he , ef f _he , dev_WNl , 
dev_ef f1) ; 

37 printf (”\n\nWith intercooling \t\t\t\t %f \t\t\t %f 

\t\t\t\t\t %f \t\t\t %f” ,WN_ic,eff_ic,dev_WN2 , 
dev_ef f2 ) ; 

38 printf (”\n\nWith Heat Exchanger & Intercooling \t\t 

%f \t\t\t %f\t\t\t\t\t %f \t\t\t %f”,WN_iche, 
eff_iche,dev_WN3,dev_eff3); 


Scilab code Exa 5.13 Comparison of Carnot efhciency with Brayton efh- 
ciency 

1 clc ; 
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2 Tl=27 + 273; // Temperature at state 1 in kelvin 

3 T3 = 827 + 273; // Temperature at state 3 in kelvin 

4 Cp = 1.005; // Specitic lieat at constant pressure in 

kJ/kg K 

5 r = 1.4; // Specific lieat ratio 

6 

7 t = T3/Tl; // Temperature ratio 

8 Wmax=Cp*((T3*(1-1/ sqrt (t)))-T1*( sqrt (t)-1)); // 

Maximum work 

9 ef f _wmax = (1 -1/ sqrt (t) ) ; // Efficiency of brayton 

cycle 

10 Tmax=T3; Tmin=Tl; 

11 ef f _carnot = (Tmax-Tmin) /Tmax ; // Carnot efficiency 

12 disp (”kJ/kg of air ” , Wmax , ”Maximum net work per kg 

o f a i r = ” ) ; 

13 disp (”%” , ef f _wmax*100 , ” Brayton cycle efficiency = 

) ; 

14 disp (”%” , ef f .carnot *100 , ” Carnot cycle efficiency = 

”) ; 


Scilab code Exa 5.15 Calculation of Improvement in Efficiency 

1 clc ; 

2 p1=1; // Pressure at state 1 in bar 

3 Tl=300; // Temperature at state 1 in kelvin 

4 p4=5; // Pressure at state 4 in bar 

5 T5 = 1250; // Temperature at state 5 in kelvin 

6 Cp = 1.005; // Specific heat at constant pressure in 

kJ/kg K 

7 r = 1.4; // Specific heat ratio 

8 

9 rp=p4/pl; // pressure ratio 

10 p2 = sqrt (rp); // Because of perfect intercooling 

11 c 1 = p2 ~ ((r -1 ) /r) ; 

12 T2 = Tl*cl; // Temperature at state 2 
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T4=T2; T3=T1; 


Wc 1 = Cp* (T2-T1) ; // Work of compressor 1 
Wc=2*Wcl ; // net work of compressor 

WTl=Wc; 

T6 = T5 - (WTl/Cp) ; // Temperature at state 6 
p5_p6=(T5/T6)“(r/(r-l)); // Pressure ratio 
p6=rp/p5_p6; // Pressure at state 6 

p7=pl; T7=T5;p8=p6; 

T8 = T7* (p7/p8) ~ ((r-1)/r) ; // Temperature in state 8 

WT2 = Cp*(T7-T8) ; // Turbine 2 work 

q=Cp* (T5-T4)+Cp* (T7-T6) ; // Heat supplied 

eff = WT2/q; // Efficiency of the cycle 

// With regenerator 

T9=T8; 

q_withregen = Cp * ((T5-T9) + (T7-T6 ) ) ; // Heat supplied 
with regenerator 

eff_withregen = WT2/q_withregen; // Efficiency of the 
cycle with regenerator 
I_eff=(eff_withregen-eff)/eff_withregen; // 

Percentage improvement in efficiency 

disp (”%” , ef f *100 , ” Efficiency of the cycle = ”,”kJ/ 
kg” , q , ” Heat supplied = ” , ” kJ/kg” , WT2 , ” Work of 
turbine = ”,”(i). Without regenerator ”); 
disp (”%” , ef f _withregen* 100 , ” Efficiency of the cycle 
= ”,”kJ/kg (roundoff er ror ) ” , q_withregen , ” Heat 
supplied = ”,”(ii). With regenerator” ); 

disp (”%” , I_ef f *100 , ” Percentage improvement in 
efficiency = ”); 


Scilab code Exa 5.16 Calculation of Efficiency ratio of the power plants 
clc ; 
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p1=1; // pressure at inlet in bar 

Tl=27 + 273; // Temperature at inlet in kelvin 

T4 = 1200; // Maximum temperature iu ke 1 vin 

t = T4/Tl; // Temperature ratio 

r = 1.4; // Specific lieat ratio 

rp = t; 

c = rp~((r-1)/r) ; 

x=(1 - sqrt (c)/rp)/(1-c/rp); 

ef f 2_1=x; 

rl = sqrt (rp) ; 

r2 = rl; r3 = rl ; r4 = rl; 

disp ( ef f 2_1 , ” Ef f i c iency ratio of power plants = ”); 
disp (r4,”pressure ratio of LPT = ” , r3 , ” pressure 

ratio of HPT = ”, r2 ,”pressure ratio of HPC = ”,rl 
,”pressure ratio of LPC = ”); 


Scilab code Exa 5.19 Determination of Net power output 
clc ; 

ra = 30; // Mass flow rate in kg/s 

p 1 = 1 ; // pressure of air at compressor inlet in bar 
Tl=273 + 15; // Temperature of air at compressor inlet 

i n k e 1 v i n 

p2 = 10.5; // Pressure of air at compressor outlet 
T_R = 420; // Temperature rise due to combustion in 

k e 1 v i n 

p4 = 1.2; // Pressure at turbine outlet in bar 

Cp = 1.005; // Specific heat at constant pressure in 

kJ/kg K 

r = 1.4; // Specific heat ratio 

T2 = T1 * (p2/pl) ~ ((r-1)/r) ; // Temperature at state 2 
T3 = T2 + T_R; // Temperature at state 3 
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13 p3 = p2; 

14 T4=T3/(p3/p4)~((r-l)/r); 

15 Wc = m*Cp* (T2-T1) ; // Compressor work 

16 WT = m*Cp* (T3-T4) ; // Turbine work 

17 WN = WT-Wc; // Net work output 

18 Q=m*Cp* (T3-T2) ; // Heat supplied 

19 eff_th = WN/Q; // Thermal efficiency 

20 

21 disp (”%” , ef f _th*100 , ” Thermal efficiency = ” , ”kW ( 
roundoff er r or ) ” , WN , ” Power output = ” , ”kW” , Q , ” 
Heat supplied = ”); 
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Chapter 6 

Practical Cycles and their 
Analysis 


Scilab code Exa 6.1 Calculation of Net power and overall efficiency of in- 
stallation 

clc ; 

p01=1; // Pressure at state 1 in bar 
T01=30 + 273; // Temperature at state 1 in kelvin 

p02 = 6; // Pressure of air after compressed in bar 
eff_c=0.87; // Isentropic efficiency of compressor 

T03=700 + 273; // Temperature at state 3 in kelvin 

eff_T = 0.85; // Isentropic efficiency of the turbine 

CV=43.1; // calorific value of fuel in MJ/kg 
ma = 80; // Mass flow rate of air in kg/rnin 

Cpa = 1.005;// Specific lieat of air at constant 
pressure in kJ/kg K 

Cpg = 1.147;// Specific heat of fuel at constant 
pressure in kJ/kg K 
rg = 1.33;// Specific heat ratio of fuel 
r = 1.4; // Specific heat ratio of air 

T_02 = T0 1 * ( p02/pO 1 ) ~ (( r-1 ) /r ) ; // from T—S diagram 
T02 = T01 + (T_02-T01) /eff _c ; // Temperature after 
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compression 

17 // Neglecting the addition of fuel in the combustion 

chamber we have mf+ma=ma 

18 mf=(ma/60) *Cpg*(T03-T02)/(CV*10~3) ; 

19 ma_mf = (ma/60) * (l/mf ) ; // Air fuel ratio 

20 A_F = ma_mf ; 

21 p04 = p01;p03 = p02 ; 

22 T_04 = T03*(p04/p03)“((rg-1)/rg) ; 

23 T04=T03-eff_T*(T03-T_04); 

24 WN=(ma/60)*Cpg*(T03-T04)-(ma/60)*Cpa*(T02-T01 ); // 

The net power of installation 

25 ef f _th = WN / (mf * C¥ * 10 ~ 3) ; // The overall thermal 

efficiency 

26 

27 disp (A_F,”(i) . Air fuel ratio of the turbine gases = 

”) ; 

28 disp (”K” , T04 , ” ( i i ) . The final temperature of exhaust 

gases = ”); 

29 disp (”kW’ , WN , ” ( i i i ) . The net power of installation = 

”) ; 

30 disp (”%” , ef f _th*100 , ” ( iv ) . The overall thermal 

efficiency = ”); 


Scilab code Exa 6.2 Calculation of ratio of compressor to turbine work 

1 clc ; 

2 p 01 = 1; // Air inlet pressure in bar 

3 T01=7 + 273;// Air inlet temperature in kelvin 

4 p02=4; // Pressure at state 2 in bar 

5 ef f _c =0.82 ; // Isentropic efficiency of compressor 

6 T03=800 + 273; // Maximum temperature at the turbine 

i n 1 e t i n k e 1 v i n 

7 eff_T = 0.85; // Isentropic efficiency of the turbine 

8 CV=43.1; // calorific value of fuel in MJ/kg 

9 Cpa = 1.005;// Specific lieat of air at constant 
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pressure in kJ/kg K 

Cpg = 1.147;// Specific lieat of fuel at constant 
pressure in kJ/kg K 
rg = 1.33;// Specific heat ratio of fuel 
r = 1.4; // Specific lieat ratio of air 

LS =0.85; 

mf = l; // Let assume rnass of fuel to be 1 kg 

T_02 = T01 * (p02/p01 ) “ ( (r-1)/r) ; // from T—S diagram 
T02 = T01 + (T_02-T01 ) /ef f _c ; // Temperature after 
compression 

Wc = Cpa* (T02-T01 ) ; // Work of compression 
Q = Cpg * ( T03-T02 ) ; // Heat supplied 
p04=p01;p03=p02; 

T_04 = T03*(p04/p03)“((rg-l)/rg) ; 
T04=T03-eff_T*(T03-T_04); 

WT = Cpg* (T03-T04) ; // Turbine work 
WN = WT-Wc; // Net work done 

ef f _th = WN/( Q/LS ) ; // The thermal efficiency 
ma_mf = (LS*CV* 10“3/Q)-1; // AIR FUEL ratio 
ma = mf *ma_mf; 

sf c = (3600/ ( ma_mf *WN) ) ; // specific fuel consumption 
Wc_WT= ( Wc*ma) /( WT *( ma+mf )); // work ratio 

disp (”kJ/kg of air ” , Wc , ” ( i ) . Compressor work = ”); 

disp (”kJ/kg of air ” , Q , ” ( i i ) . Heat supplied = ”); 

disp (” k J/ kg o f air” , WT , ” ( iii ) . Turbine work = ” ) ; 

disp (” kJ /kg o f air” , WN ,” ( iv) . Net work = ”); 

disp (”%” , ef f _th*100 , ” (v) . Thermal Efficiency =”); 
disp (ma_mf , ” ( vi) . Air / Fuel ratio = ”) 
disp (” kg/kWh” , sf c , ” ( vii ) . Specific fuel consumption 
=”) ; 

disp (Wc_WT , ” ( v i i i ) . Ratio of compressor work to 
turbine work = ”); 
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Scilab code Exa 6.3 Calculation of effect of pressure loss 
clc ; 

eff_c=0.82; // Isentropic efficency of the 

compressor 

eff_T = 0.85; // Isentropic efficency of the turbine 

eff_m = 0.99; // Mechanical transmission efficiency 

r p = 7 ; // Pressure ratio 

T03 = 1000; // Maximum cycle temperature in kelvin 

ef f _comb=0.97 ; // Combustion efficiency 

CV=43.1; // Calorific value in MJ/kg 
ma = 20; // Air mass flow rate in kg/s 

eff_reg=0.75; // Regenerator effectiveness 

del_P=0.1; // Regenerator gas side pressure loss in 

bar 

T01=327; // Ambient temperature in kelvin 

p01 = 1; // Ambient pressure in bar 

Cpa = 1.005;// Specific heat of air at constant 
pressure in kJ/kg K 

Cpg = 1.147;// Specific heat of fuel at constant 
pressure in kJ/kg K 
rg = 1.33;// Specific heat ratio of fucl 
r = 1.4; // Specific lieat ratio of air 

//( i ) . With Regeneration and pressure loss 

T_02 = T01*(rp)~((r-l)/r) ; 

T02 = T01 + (T_02-T01)/eff_c ; 
p04=p01+del_P; 
p03=rp/pO1; 

T_04 = T03*(p04/p03)~((rg-l)/rg) ; 
T04_l=T03-eff_T*(T03-T_04); 

T05=T02+eff_reg*(T04_l-T02); 

mf1 = (ma*Cpg*(T03-T05))/(CV* 10~3*eff_comb) ; // By 

neglecting the effect of change in mass flow rate 
due to mf in combustion chamber 

p03_p04_l=p03/p04; 

WTl = (ma+mf1 )* Cpg *( T03-T04_1 ); // Turbine work 
WNl=(ma+mfl)*Cpg*(T03-T04_l)-(ma*Cpa*(T02-T01)/eff_m 
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) ; // Net work output 

sf c 1 =mf 1 *3600/WN1 ; // Specifc fuel consumptiou 
ef f _th 1 =WN1 / (mf 1 *CV* 10~3) ; // Thermal efficiency 


//( i i ) . Without Regenerator and without pressure loss 

p04=p01; 

T_04 = T03*(p04/p03)“((rg-l)/rg) ; 
T04_2=T03-eff_T*(T03-T_04); 

mf2=(ma*Cpg*(T03-T02))/(CV*10~3*eff_comb); 
WT2=(ma*Cpg*(T03-T04_2)); 

WN2=(ma*Cpg*(T03-T04_2))-(ma*Cpa*(T02-T01)/eff_m); 

// Net work output 
p03_p04_2=p03/p04; 

sfc2 = mf 2*3600/WN 2; // Specific fuel consumption 
ef f _th2 = WN2/ (mf 2*CV* 10~3) ; // Thermal efficiency 


//( i i i ) . With Regenerator and without pressure loss 

T_02 = T01*(rp)~((r-l)/r) ; 

T02 = T01 + (T_02-T01)/eff_c ; 
p04=p01; 
p03=rp/pO1; 

T_04 = T03*(p04/p03)~((rg-l)/rg) ; 
T04_3=T03-eff_T*(T03-T_04); 

T05 = T02 + eff_reg*(T04_3-T02) ; 

WT3=(ma*Cpg*(T03-T05)); 

mf3 = (ma*Cpg *( T03-T05 ))/( CV * 10 ~ 3*eff_comb ); // By 

neglecting the effect of change in mass flow rate 
due to mf in combustion chamber 

p03_p04_3=p03/p04; 

WN3=(ma+mf3)*Cpg*(T03-T04_3)-(ma*Cpa*(T02-T01)/eff_m 

) ; // Net work output 

sf c3 = mf 3 *3600/WN3 ; // Specifc fuel consumption 
ef f _th3 = WN3/ (mf 3*CV* 10~3) ; // Thermal efficiency 


26 


64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

1 

2 

3 

4 


printf (” Quantities \t\t\t \ t \ t Regener at or \t\t\t\t\t 
\ t Without” ) ; 

printf (” \n\t\t\t\twith Del_P \ t \ twit hout Del_P\t\t\t 
\ tregenerator and Del_P”); 

printf (”\n\t\t\t\t ( roundoff error )\t ( roundoff error 
)\t\t\t(roundoff error)”); 
printf (”\n\n P03/P04\t\t\t%f\t\t%f\t\t\t\t\t%f” , 
p03_p04_l,p03_p04_3,p03_p04_2); 
printf (”\n\nT04 (K)\t\t\t \t%f\t\t%f\t \t\t\t\t%f” , 
T04_ 1 , T04_3,T04_2 ) ; 

printf (”\n\nmf (kg/s) \ t\ t\ t%f\t\t%f\t\t\t\t\t%f” , 

mf 1 , mf 3 , mf 2) ; 

printf (”\n\nWT (kW)\t\t\t \t%f\t\t%f\t \t\t\t\t%f” , 

WTl , WT3 , WT2 ) ; 

printf (”\n\nsfc (kg/kW h) \t\t\t%f\t\t%f\t\t\t\t\t%f 

”,sfcl,sfc3,sfc2); 

printf (”\n\nefficiency (in percentage ) \t%f\t\t%f\t \ 
t\t\t\t%f”,eff_th1*100,eff_th3*100,eff_th2*100); 

printf (”\n\nAs can be seen from the tablc that 
pressure loss plays a major role in the 
efficiency tlian the regenerator . \n\nHence , more 
care should be taken in the design to have 
minimum pressure 1 oss . ” ) ; 


Scilab code Exa 6.4 Calculation of net power out SFC and overall Efh- 
ciency 

clc ; 

eff_c=0.8; // Isentropic efficiency of compression 

each stage 

eff_CT=0.88; // Isentropic efficiency of compressor 

t u r b i n e 

eff_PT=0.88; // Isentropic efficiency of power 


27 



5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 


t u r b i n e 

ef f _trans =0.98 ; // Turbine to compressor 

transmission efficiency 

rp = 3; // Pressure ratio in eacli stage of compression 

T08 = 297; // Temperature after intercooler in kelvin 

ma = 15; // Air mass flow in kg/s 
eff_reg = 0.8; // Regenerator effect i veness 

del_P=0.1; // Regenerator gas side pressure loss in 

bar 

T01=327; // Ambient temperature in kelvin 

p01 = 1; // Ambient pressure in bar 

T03 = 1000; // Maximum cycle temperature in kelvin 

CV = 43.1; // Calorific value in MJ/kg 
Cpa = 1.005;// Specific lieat of air at constant 

pressure in kJ/kg K 

Cpg = 1.147;// Specific lieat of fuel at constant 
pressure in kJ/kg K 
rg = 1.33;// Specific heat ratio of fuel 
r = 1.4; // Specific heat ratio of air 

p03=rp~2; // Pressre at state 3 in bar 

T_07=T01*(rp)~((r-l)/r); 

T07 = T01 + (T_07-T01)/eff_c ; 

WLPC = ma*Cpa* (T07-T01 ) ; // Work of low pressue 
compressor 

T_02 = T08*(rp)~((r-l)/r) ; 

T02 = T08 + (T_02-T08)/eff_c ; 

WHPC=ma*Cpa*(T02-T08); 

WC = WLPC + WHPC ; // Compressor work 
WCa = WC/ef f _trans ; // Actual compressor work 
// Neglecting effect of mf 
T09=T03-(WCa/(ma*Cpg)); 

T_09=T03-(T03-T09)/eff_PT; 
p09 = p03/(T03/T.09)"(rg/(rg-1)) ; 
p04=p01+del_P; 

T_04 = T09*(p04/p09)“((rg-1)/rg) ; 

T04=T09-eff_PT*(T09-T_04); 

WTP = ma*Cpg* (T09-T04) ; // Work output of power 
t u r b i n e 
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T05=T02+eff_reg*(T04-T02); 

mf = (ma*Cpg*(T03-T05))/(CV*10~3) ; 

sf c=mf *3600/(WTP) ; //Specifc fuel consumption 

ef f _th = WTP/ (mf *CV* 10~3) ; // Thermal efficiency 


disp (”kW (roundoff err or ) ” , WTP , ” Work output of 
power turbine = ”); 

disp (”kg/kW h” , sf c , ” Specifc fuel consumption = ”); 
disp (”%” , ef f _th*100 , ” Thermal efficiency =”); 


Scilab code Exa 6.5 Calculation of the thermal efficiency and air rate 
clc ; 

Wplant = 1850; // Plant work output in KW 

p01 = 1; // Ambient pressure in bar 

T01=27 + 273; // Ambient temperature in kelvin 

T03=720 + 273; // Maximum cyc 1 e temperature in ke 1 vin 

rp=2.5; // Pressure ratio 

eff_T = 0.80; // Turbine and compressor efficiency 

eff_reg=0.75; // Regenerator effectiveness 

ef f _comb=0.98 ; // Combustion efficiency 

CV=43.1; // Calorific value in MJ/kg 
del_p=0.03; // Pressure drop 

p02=6.25; // Pressure in bar 

Cpa = 1.005;// Specific lieat of air at constant 
pressure in kJ/kg K 

Cpg = 1.147;// Specific heat of fuel at constant 
pressure in kJ/kg K 
rg = 1.33;// Specific lieat ratio of fucl 
r = 1.4; // Specific heat ratio of air 

T_07 = T01*rp~((r-l)/r) ; 

T07 = T01 + (T_07-T01)/eff_T ; 

T02 = T07; 
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WLPC = Cpa* (T07-TOl) ; // Work of low pressure 
compressor 

WHPT = WLPC ; 

T09=T03-WHPT/Cpg; 

T_09 = T03-(T03-T09)/eff_T ; 
p03=(l-del_p)~2*p02 
p09 = p03/(T03/T_09)~(rg/(rg-1)) ; 
plO=p09*(1-del_p); 

T10=T03; 
p04=p01+del_p; 

T_04=T10*(p04/pl0)~((rg-l)/rg); 
T04=T10-eff_T*(T10-T_04); 

Wlpt=Cpg*(T10-T04); 

WN=Wlpt-WHPT; 
ma=Wplant/WN; 

T05=T02+eff_reg*(T04-T02); 

Q = Cpg*(T03-T05 + T10-T09) ; 
eff_th=WN/Q; 

WHPT_l=ma*WHPT; 

Wlpt_l=ma*Wlpt; 

mf=ma*Q*3600/(eff_comb*CV*10~3) ; 
sfc=mf/Wplant; 

disp (”K”,T_07,” T_07 = ”); 
disp (”K” ,T07, ” T07 = ”); 
disp (”K”,T09, ” T09 = ”); 
disp (”K” ,T_09,” T_09 = ”); 
disp (”K” , T_04,” T_04 = ”); 
disp (”K”,T04, ” T04 = ”) ; 
disp (”K” ,T05, ” T05 = ”) ; 
disp (” bar” ,p03 ,”P03 = ”) ; 
disp (” bar” ,p09 ,”P09 = ”) ; 
disp (” bar” ,pl0 ,”P10 = ”) ; 
disp (” kg/s ” , ma , ” Mass flow rate = ”); 
disp (”%” , ef f _th*100 , ” The overall efficiency =”); 
disp (” kg o f f u e 1 /kW h” ,sfc ,” Specific fuel 
consumption = ” ) ; 
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Scilab code Exa 6.6 Calculation of Compressor efficiency and the temper- 
ature ratio 

clc ; 

rp=11.3137; // Pressure ratio 

WN = 0; // Net work output 

Q=476.354; // Heat added per kg of air mass in kJ 

T01=300; // Inlet air total temperature in k e 1 v i n 

eff_T = 0.71; // turbine efficiency 

Cpa = 1.005;// Specific heat of air at constant 

pressure in kJ/kg K 

Cpg = 1.147;// Specific lieat of fuel at constant 
pressure in kJ/kg K 
rg = 1.33;// Specific heat ratio of fuel 
r = 1.4; // Specific lieat ratio of air 

T_02 = T01*rp"((r-l)/r) ; 

T03_T02=Q/Cpa; 

T03_T_04 = rp~((r-1)/r) ; 

T04_T03=l-(eff_T*(l/T03_T_04)*(T03_T_04-l)); 
T04=T01+(T03_T02); 

T03=T04/T04_T03; 
t = T03/T01; // Temperature ratio 
T02=T03-T03_T02; 

ef f _C = (T_02-T01) / (T02-TO1) ; // Compressor efficiency 

disp (”%” , ef f _C*100 , ” Compressor Efficiency = ”,); 
disp (t , ” Temperature ratio = ”); 


Scilab code Exa 6.7 Calculation of suitable pressure ratio 
clc ; 
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eff_C = 0.7042; // Efficiency of the compressor 

eff_T=0.71; // Efficiency of the turbine 
Q=476.354; // Head added in kJ/kg 

WR = 0.0544; // Work ratio 

T01=300;// Total inlet temperature in kelvin 
Cpa = 1.005;// Specific lieat of air at constant 
pressure in kJ/kg K 

Cpg = 1.147;// Specific heat of fuel at constant 
pressure in kJ/kg K 
rg = 1.33;// Specific lieat ratio of fucl 
r = 1.4; // Specific heat ratio of air 

c_t = (l-WR)*(eff_T*eff_C) ; 

t = ((Q/(Cpg*T01))+l-l/eff_C)/(1-c_t/eff_C) ; // 
Temperature ratio 
c = c_t *t ; 

rp=c“(r/(r-l)); // Pressure ratio 

disp (rp,”Pressure ratio = ”); 
disp (t , ” Temperature ratio = ”); 


Scilab code Exa 6.8 Calculation of minimum temperature ratio 
clc ; 

WR = 0.3; // Work ratio 

rp=12; // Pressure ratio 

t=4; // Temperature ratio 

Cpa = 1.005;// Specific heat of air at constant 
pressure in kJ/kg K 

Cpg = 1.147;// Specific heat of fuel at constant 
pressure in kJ/kg K 
rg = 1.33;// Specific lieat ratio of fuel 
r = 1.4; // Specific lieat ratio of air 

c=rp~((r-1)/r); 


32 



11 eff_C_T=l/((l-WR)*t/c); 

12 tmin=c/eff_C_T; 

13 ef f =1 -1/c ; 

14 

15 disp (tmin , ” Minimum Temperature ratio = ”); 

16 disp (”%” , ef f *100 , ” Efficiency = ”); 


Scilab code Exa 6.9 Calculation of Isentropic efficiency of Turbine 

1 clc ; 

2 eff_pe=0.85; // Polytropic efficiency of the 

compressor 

3 T_02_T01=2; 

4 Cpa = 1.005;// Specific lieat of air at constant 

pressure in kJ/kg K 

5 Cpg = 1.147;// Specific lieat of fuel at constant 

pressure in kJ/kg K 

6 rg = 1.33;// Specific heat ratio of fuel 

7 r = 1.4; // Specific heat ratio of air 

8 

9 rc = (T_02_T01)~(r/(r-1)) ; 

10 eff_C = (T_02_T01-1)/(((rc~(((r-1)/r)*(1/eff_pe))) -1)) 

; // Compressor efficiency 

11 eff_T=(l-(l/rc)“(eff_pe*(r-l)/r))/(l-(l/rc)~((r-l)/r 

)); // Turbine efficiency 

12 

13 

14 disp (”%” , ef f _C*100 , ” Isentropic compressor 

efficiency = ”); 

15 disp (”%” , ef f _T*100 , ” Isentropic Turbine efficiency 

= ”); 
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Scilab code Exa 6.10 Plotting yariation of Isentropic efficiency over a range 
of pressure ratio 

1 clc ; 

2 eff_C = 0.85; // Isentropic efficiency of the 

compressor 

3 rp=4; // Pressure ratio 

4 r = 1.4; // specific heat ratio 

5 eff_pc = (((r-1) /r)*log (rp))/log (((rp~((r-1)/r)-1)/ 

eff_C)+l); 

6 disp (”%” , eff _pc*100 , ” Polytropic efficiency = ”); 

7 disp (”variation of compressor efficiency witli 

compression ratio is shown in windowl”); 

8 xset ( ’window ’ , 1) ; 

9 function eff_c=f(rc) 

10 eff_c = (rc~0.286-l)/(rc~0.326-l) ; 

11 endfunction 

12 rc = linspace (2,10,4); 

13 plot (rc , f ) ; 

14 title (”variation of compressor efficiency with 

compression ratio” ,” fontsize” ,4,” color” , ” blue ” ) ; 

15 xlabel(”compression ratio (rc)”,”fontsize”,4,”color” 

,”blue”); 

16 ylabel (”Compressor e f f i c i en cy ” , ” f o n t si z e ” , 4 , ” color ” 

,”blue”); 


Scilab code Exa 6.11 Calculation of the power output thermal efficiency 
and the heat exchanger effectiveness 

1 clc ; 

2 eff_pe=0.88; // Compressor and turbine polytropic 

efficiencies 

3 T01=310; // Temperature at LP compressor inlet in 

k e 1 v i n 

4 p01 = 14; // Pressure at LP compressor inlet in bar 
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rp = 2; // Compressor pressure ratio 

T03 = 300;// Temperature at HP compressor inlet in 
k e 1 v i n 

m = 180; // Mass flow of Helium in kg/s 

Q = 500; // Heat inpnt to gas tnrbine in MW 

T07=700; // Helium Temperature at entry to reactor 

c h a n n e 1 s i n k e 1 v i n 

P_precoller=0.34; // Pressure loss in pre—cooler and 

intercooler in bar 

P_loss_HE=0.27; // Pressure loss in heat exchanger 
in bar 

P_loss_RC = 1.03; // Pressure loss in reactor channel 
in bar 

eff_pc=0.88; // Polytropiic efficiency 

Cp =5.19; // Specific heat at constant pressure in 
kJ/kg K 

r = 1.66; // Specific heat ratio 

n_l_n=((r-l)/r)*(1/eff_pc) ; 

T02=T01*rp"n_l_n; 

T04=T03*rp~n_l_n; 

T05=((Q*10~3)/(m*Cp))+T07; 

T_press_loss=P_precoller+P_loss_HE+P_loss_RC; // 

Total pressure loss 
p05=56-T_press_loss; 
p06=p01+P_precoller+P_loss_HE; 
n__l_n = eff_pc*((r-l)/r) ; 

T06=T05/(p05/p06)~n__l_n; 

WC = m*Cp* ((T02-T01) + (T04-T03) ) ; // Work of compressor 

WT = m*Cp* (T05-T06 ) ; // Work of Turbine 

WN = WT-WC; // Net work output 

ef f _th = WN/ ( Q* 10“3) ; // Efficiency 

eff=(T07-T04)/(T06-T04); // Effectiveness 

disp (”MW (roundoff err or ) ” , WN/1000 , ” Power output 
= ”); 

disp (”% (roundoff error ) ” , ef f _th*100 , ” Thermal 
efficiency = ”); 
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disp (”% (roundoff error ) ” , ef f *100 , ” Effectiveness 
= ”) ; 


Scilab code Exa 6.12 Calculation of the isentropic efficiency of the turbine 
and the mass flow of air 

clc ; 

rp=4; // Pressure ratio 

WN = 1500; // Net work output in kW 

T01=25 + 273; // Inlet temperature iu k e 1 v i n 

p 01 = 1; // Inlet pressure in bar 

p03=4; // Turbine inlet pressure in bar 

T03=700 + 273;// turbine inlet temperature in kelvin 
eff_c=0.85; // Compressor efficiency 

eff_over=0.21; // Overall efficiency 
Cp = 1.005;// Specific heat of air at constant 
pressure in kJ/kg K 
r = 1.4; // Specific heat ratio of air 

T02=T01+T01*(rp"((r-1)/r)-1)/eff_c; 

Q=WN/eff_over; 
m=Q/(Cp*(T03-T02)); 

Wn = WN/m; // Net work per kg 
T04=T03-T02+T01-(Wn/Cp); 

T_04=T03/rp~((r-1)/r); 
eff_T=(T03-T04)/(T03-T_04); 

disp (” kg/ s ” , m , ” Mass flow rate = ”); 
disp (”%” , ef f _T*100 , ” Isentropic efficiency of the 
Turbine = ” ) ; 


Scilab code Exa 6.13 Determination of pressure of the gas entering the 
low pressure turbine 
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clc ; 

rp=4; // Pressure ratio 

eff_c=0.86; // Compressor efficiency 

ef f _Thp =0.84 ; // High pressure turbine efficiency 

eff_Tlp=0.8; // Low pressure turbine efficiency 

eff_M = 0.92; // Mechanical efficiency 

T03 = 660 + 273; // in kelvin 

T05 = 625 + 273; // I n k e1vin 

T01 = 15 + 273; // Inlet temperature in k e 1 v i n 

p 01 = 1; // Inlet pressure in bar 

Cp = 1.005; // Specific lieat of air at constant 
pressure in kJ/kg K 
r = 1.4; // Specific lieat ratio of air 

eff= 0.75; // Heat exchanger effect i veness 

T_02=T01*(rp)~((r-l)/r); 

T02=((T_02-T01)/eff_c)+T01; 

T04 = T03-((T02-T01)/eff_M) ; 

/ / I n HP t u r b i n e 

T_04 = T03-((T03-T04)/ef f _Thp); 

p_04 = rp/(T03/T_04)~(r/(r-l)) ; 

II In LP t u r bin e 
p05=p_04;p_06=p01; 

T_06=T05/(p05/p_06)~((r-l)/r); 
T06=T05-(eff_Tlp*(T05-T_06)); 

T07=T02+eff*(T06-T02); 

Q=Cp*(T03-T07+T05-T04); 

Wc=Cp*(T02-T01); 

WT=Cp*(T03-T04+T05-T06); 
eff_th=(WT-Wc)/Q; 

disp (” bar” , p_04 , ” ( i ) . Pressure of gas entering low 
pressure turbine = ”); 

disp (”%” , ef f _th*100 , ” Overall efficiency =”); 
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Scilab code Exa 6.14 Calculation of pressure ratio and cycle efficiency 


clc ; 

T01=38 + 273; // Inlet temperature of compressor in 

k e 1 v i n 

eff_c=0.82; // Compressor efficiency 

T03 = 650 + 273; // Turbine inlet temperatnre in kelvin 

eff_T = 0.8; // Turbine efficiency 

Cpa = 1.005;// Specific heat of air at constant 
pressure in kJ/kg K 

Cpg = 1.147;// Specific heat of fuel at constant 
pressure in kJ/kg K 
rg = 1.33;// Specific heat ratio of fuel 
r = 1.4; // Specific lieat ratio of air 

t=T03/T01; 

// For maximun specific work we know that 

ropt = (sqrt (t*eff_c*eff_T))~(r/(r-l)) ; 

T_02 = T01 *ropt~((r-1)/r) ; 

T02 = T01 + (T_02-T01)/eff_c ; 

T_04=T03/ropt~((rg-l)/rg); 

T04=T03-eff_T*(T03-T_04); 

eff_th=((Cpg*(T03-T04))-(Cpa*(T02-TO1)))/(Cpg*(T03- 
T02 ) ) ; 

disp (ropt , ” Optimum pressure ratio = ”); 

disp (”%” , ef f _th*100 , ” Overall efficiency = ”); 


Scilab code Exa 6.15 Estimation of the temperature of the gases at entry 
to the turbine 

clc ; 

p 0 1 = 1 ; // Stagnation pressure at entry in bar 

pa=0.93; // Static pressure at entry in bar 

Tl = 10 + 273; // Static temperature in entry in ke 1 vin 
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p 0 2 = 6 ; // Pressure at state 2 in bar 

T02=230 + 273; // Temperature at state 2 in kelvin 

P = 5100; // Turbine output power in kW 

A = 0.1; // Compressor entry area in m' 2 

Cpa = 1.005;// Specific lieat of air at constant 
pressure in kJ/kg K 

Cpg = 1.147;// Specific heat of fuel at constant 
pressure in kJ/kg K 
rg = 1.33;// Specific lieat ratio of fucl 
r = 1.4; // Specific heat ratio of air 

R = 287; // Characteristic constant in J/kg K 
T04=460 + 273; // Exhaust pipe temperature in kelvin 

M=sqrt (((p01/pa)“((r-1)/r)-1)/((r-1)/2)) ; 
T01=Tl*(l+(r-l)/2*M“2); 

T_02=T01*(p02/p01)“((r-l)/r); 
eff_c=(T_02-T01)/(T02-T01); 
row_s=(pa*10~5)/(R*T1); 
a=sqrt (r*R*Tl); 

V=M*a; 

ra = row_s *A*V; 

T03=(P/(m*Cpg))+T04; 

disp (”%” , ef f _c *100 , ” Compressor efficiency = ”); 
disp (” kg/s ” , m , ” Mass flow rate = ”); 
disp (”K (roundoff error ) ” , T03 , ” Turbine inlet 
stagnation temperature = ”); 


Scilab code Exa 6.16 Calculation of efficiency and work ratio of modern 
turbines and old turbines 

clc ; 

T01=27 + 273; // Inlet temperature in kelvin 
p 01 = 1; // Inlet pressure in bar 

rp=3; // Pressure ratio 
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Cpa = l . 005; // Specific heat of air at constant 
pressnre in kJ/kg K 

Cpg = 1.147;// Specific lieat of fuel at constant 
pressure in kJ/kg K 
rg = 1.33;// Specific heat ratio of fuel 
r = 1.4; // Specific lieat ratio of air 

R = 287; // Characteristic constant in J/kg K 

T_02 = T01*rp"((r-l)/r) ; 

// Turbines 70 years ago 

eff_c=0.65; // Compressor efficiency 

eff_T = 0.7; // Turbine efficiency 

T03 =700 + 273; // in kelvin 

T02=T01*(l+((rp“((r-l)/r)-l)/eff_c)); 

T04=T03*(l-eff_T*(l-(l/rp“((rg-l)/rg)))); 

eff_th=(Cpg*(T03-T04)-Cpa*(T02-T01))/(Cpg*(T03-T02)) 

J 

WR=(Cpg*(T03-T04)-Cpa*(T02-T01))/(Cpg*(T03-T04)); 

disp (WR , ” Work ratio = ” , ef f _th*100 , ” The Efficiency 
= ” ,” ( i ) .70 years ago”); 

//Modern turbines 

eff_c=0.85; // Compressor efficiency 

eff_T = 0.9; // Turbine efficiency 

T03=1000+273; // in kelvin 
T02 = T01 + (T_02-T01)/eff_c ; 

T_04=T03/rp"((rg-1)/rg); 

T04=T03-eff_T*(T03-T_04); 

Wc=Cpa*(T02-T01); 

WT=Cpg*(T03-T04); 

WN=WT-Wc; 

eff_th=WN/(Cpg*(T03-T02)); 

WR=WN/WT; 

disp (WR , ” Work ratio = ef f _th * 100 ,” The 

Efficiency = , ” ( i i ) . Modern turbines”); 
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Scilab code Exa 6.17 Determination of necessary mass flow rate 
clc ; 

r p = 7; // Pressure ratio 

T03 = 1000; // Maximum temperature in ke 1 vin 

eff_c=0.85; // Compressor efficiency 

eff_T = 0.9; // Turbine efficiency 

T01=288; // Air entering temperature in k e 1 v i n 

PN = 750; // Power output in kW 

Cpa = 1.005;// Specific lieat of air at constant 
pressure in kJ/kg K 

Cpg = 1.147;// Specific lieat of fuel at constant 
pressure in kJ/kg K 
rg = 1.33;// Specific heat ratio of fuel 
r = 1.4; // Specific heat ratio of air 

R = 287; // Characteristic constant in J/kg K 

// Actual cycle 

T02=T01*(l+((rp“((r-l)/r)-l)/eff_c)); 

T04=T03*(l-(eff_T*(l-(l/rp~((r-1)/r))))); 
WN_a=(Cpa*(T03-T04)-Cpa*(T02-T01)); 
eff_th = WN_a/(Cpa*(T03-T02)) ; 

disp (”%” ,eff_th*100 ,”The Efficiency =”,”kJ/kg”, 
WN_a , ” Net work = ” , ” ( i ) . Actual cycles”); 

// Ideal cycle 

WN = Cpa*((T03*(l-(l/rp~((r-1)/r))))-TO1 *((rp~((r-1)/r 
)-l))); 

eff_th=l-(l/rp~((r-l)/r)); 
ma=PN/WN_a; 

disp (” kg / s ” ,ma,” Mass flow rate eff_th * 100 , 

” The Efficiency = ” , ” kJ/kg” , WN , ” Net work = ”,”(ii 
).Ideal cycles”); 
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Scilab code Exa 6.18 Estimation of intermediate pressure and tempera- 
ture between the two turbine stages 

clc ; 

ra = 5; // Mass flow rate in kg/s 

p01 = 1; // Pressure at state 1 in bar 

p02=5; // Pressure at state 2 in bar 

ef f _c =0.85 ; // Compressor efficiency 
ef f _Thp =0.87 ; // Higli pressure turbiue efficiency 

eff_Tlp=0.82; // Low pressure turbiue efficiency 

T03 = 675 + 273; // HP turbine inlet temperature in 

k e 1 v i n 

eff=0.7; // Effectiveness of the heat exchanger 

T01 = 15 + 273; // Temperature at state 1 in kel vin 
Cpa = 1.005;// Specific heat of air at constant 
pressure in kJ/kg K 
r = 1.4; // Specific heat ratio of air 

R = 287; // Characteristic constant in J/kg K 

p03=p02; 

T_02 = T01*(p02/p01)“((r-l)/r) ; 

T02 = T01 + (T_02-T01)/eff_c ; 

T04=T01-T02+T03; 

T_04 = T03-(T03-T04)/eff_Thp ; 
p04 = p03/(T03/T_04)“(r/(r-l)) ; 
p05=p01; 

T_05=T04/(p04/p05)“((r-1)/r); 

T05=T04-eff_Tlp*(T04-T_05); 

T0x=eff*(T05-T02)+T02; 

Wlpt=Cpa*(T04-T05); 

Plpt = Wlpt*m; 

Q=Cpa*(T03-T0x); 
eff_th=Wlpt/Q; 
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disp (” Intermediate pressure p04 and temperature T04 
between the two turbine stages ”); 
disp (”K” ,T04, ” To4 = ” , ” bar” ,p04,”P04 = ”); 
disp (”kW” ,Plpt, ” Power output of LP turbine = ”); 
disp (” kJ/kg” , Q , ” Heat supplied = ”); 
disp (”%” , ef f _th*100 , ” The Overall efficiency =”); 


Scilab code Exa 6.19 Determination of the percentage of the total air in- 
take that passes to the turbine and the thermal efficiency of the plant 

clc ; 

rlp=3; // Pressure ratio 

rhp = rlp ; 

eff_c=0.82; // Compressor efficiency 

T04 = 650 + 273; // Temperature at state 4 in kelvin 

T05 = 540 + 273; // Temperature at state 5 in kelvin 

eff_T = 0.87; // Efficiency of turbine 

T01 = 15 + 273; // Temperature at compressor inlet in 

k e 1 v i n 

Cpa = 1.005;// Specific heat of air at constant 
pressure in kJ/kg K 

Cpg = 1.147;// Specific heat of fuel at constant 
pressure in kJ/kg K 
rg = 1.33;// Specific heat ratio of fuel 
r = 1.4; // Specific lieat ratio of air 

T02=T01*(l+(rlp“((r-l)/r)-l)/eff_c); 
T03=T02*(l+(rhp"((r-l)/r)-1)/eff_c); 

T_06=T05/(rlp)“(2*(rg-1)/rg); 

T06=T05-eff_T*(T05-T_06); 

xl=l-((T02-T01)/(((Cpg/Cpa)*(T05-T06)-(T03-T02)))); 
x=abs (xl); 

T07=T04*(l-(eff_T*(l-(l/rhp“((rg-l)/rg))))); 
eff_th=(x*Cpg*(T04-T07))/((l-x)*Cpg*(T05-T03)+x*Cpg 
*(T04-T02)); 
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disp (”%” , (x) * 100 ,” Percentage of the total air 
intake tliat passes to the power turbine = ”); 
disp (”% ( Roundoff error ) ” , ( ef f _th) *100 , ” The 

overall efficiency = ”); 


Scilab code Exa 6.20 Calculation of the Thermal efficiency 
clc ; 

rp = 2; // Pressure ratio 

T01 = 15 + 273; // Inlet temperature in k e 1 v i n 

p 01 = 1; // Inlet pressure in bar 

T05=700 + 273; // Temperature at state 5 in kelvin 

T07=T05; 

eff_c=0.85; // compressor efficiency 

eff_T = 0.85; // Turbine efficiency 

eff=0.5; // Effectiveness of lieat exchanger 

Cp = 1.147;// Specific heat at constant pressure in kJ 
/kg K 

rg = 1.33;// Specific heat ratio of fucl 
r = 1.4; // Specific heat ratio of air 

T03=T01; 

// p02/p01=p04/p03=rp 
//p04/p01=p05/p08=rp'2 

T_02=T01*(rp)~((r-l)/r); 

T02 = T01 + (T_02-T01)/eff_c ; 

T04=T02; 

T_06=T05/rp~((rg-1)/rg); 

T06=T05-eff_T*(T05-T_06); 

T08=T06; 

T09=T04+eff*(T08-T04); 

WN=Cp*(T07-T08); 

Q=Cp*(2*T05-T06-T09); 
eff_th=WN/Q; 
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disp (” kJ/kg” ,WN,” Net work done = ”); 

disp (”%” , ef f _th*100 , ” The overall efficiency =”); 


Scilab code Exa 6.21 Calculation of cycle thermal efficiency 
clc ; 

T01=270 + 273; // Temperature at state 1 in kelvin 

T03=T01; 

p 01=1; // Inlet pressure in bar 

rp=6; // Pressure ratio 

eff_c=0.85; // Compressor efficiency 

T05 = 1150 + 273; // Temperature at inlet to expansion 

i n k e 1 v i n 

eff_T = 0.9; // Turbine efficiency 

n = 1.24 ; // Polytropic index 

R=10.05; // in kJ/kg K 

T_02=T01*rp~((n-l)/n); 

T02 = T01 + (T_02-T01)/eff_c ; 

Cv=R/(n-1); 

Cp=R+Cv; 

Wc=2*Cp*(T02-T01); 

T_06=T05/rp~((n-l)/n); 

T06=T05-eff_T*(T05-T_06); 

WT=2*Cp*(T05-T06); 

Q=Cp*(T05-T02)+Cp*(T05-T06); 

WN=WT-Wc; 
eff_th=WN/Q; 

disp (”%” , eff _th*100 , ” The Cycle efficiency =”); 


45 




1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 


Chapter 7 

Jet Propulsion Cycles and 
Their Analysis 


Scilab code Exa 7.1 Calculation of Thrust power and Heat Input 
clc ; 

CV = 43; // Calorific value of fuel in MJ/kg 

mf =0 . 18*9000/3600; // Fuel consumption in kg/s 

F = 9; // Thrust in kN 

ci=500; // Aircraft velocity in m/s 

ma = 27; // Mass of air passing through compressor in 

kg/s 

A_F = ma/mf; // Air fuel ratio 

PT = F*ci; // Thrust power 

Q = mf * ( CV * 10"3) ; // Heat supplied 

eff = PT/Q; // Overall efficiency 

disp (A_F , ” Air fuel ratio = ”); 

disp , ef f *100 , ” Overall efficiency =”); 


Scilab code Exa 7.2 Calculation of The Total Design Thrust 
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clc ; 

T03 = 1200; // Maximmn turbine inb 1 et temperature in 
k e 1 v i n 

rc=4.25; // Pressure ratio across compressor 

ma = 25; // Mass flow rate in kg/s 

eff_C = 0.87; // Isentropic efficiency of the 

compressor 

eff_T = 0.915; // Isentropic efficiency of turbine 

eff_n = 0.965; // Propclling nozzle efficiency 

eff_Tr = 0.985; // Transmission efficiency 

del_pcomb =0.21; // Combustion chamber pressure loss 

in bar 

Cpa = 1.005; // Specific heat at constant pressure of 
air in kJ/kg K 

ra = 1.4; // Specific heat ratio of air 

Cpg = 1.147; // Specific heat of fuel in kJ/kg K 

rg = 1.33; // Specific heat of fucl 

T01=293; // Ambient temperature in kelvin 

p01 = 1; // Ambient pressure in bar 

A_F=50; // Air Fuel ratio 

p02=rc/pOl; 

T02=(T01*((rc)“((ra-l)/ra)-l)/eff_C)+T01; // Actual 
temperature at state 2 
T04=T03-((Cpa*(T02-T01))/(eff_Tr*Cpg)); // 

Temperature at state 4 

rt=(l/(l-((T03-T04)/(eff_T*T03))))“(1/((rg-1)/rg)); 

// Pressure ratio across turbine 
p04=(p02-del_pcomb)/rt; // Pressure at 4 
p5=p01; 

T_5 = T04/(p04/p5) ~ ((rg-1)/rg) ; // Temperature at 5 

T5=T04-eff_n*(T04-T_5); 

c5 = sqrt (2*Cpg*10~3*(T04-T5)) ; 

F = ma*c5; // Total design thrust 

p04_pc=l/(1-((1/eff_n)*((rg-1)/(rg+1))))~(rg/(rg-1)) 
pc=p04*(l/p04_pc); 

Tc = T04*(l/p04_pc)~((rg-l)/rg) ; 

R=Cpg*10~3*(rg-l)/rg; 
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cj=sqrt (rg*R*Tc); 
row_c=(pc*10~5)/(R*Tc); 

A = ma/( row_c*c j ) ; // Area of the propelling nozzle 
d=sqrt (4*A/3.14); // Diameter of the nozzle 
pa=p01; 

Fp = (pc-pa)* 10 ~ 5*A ; // Pressure thrust 
Fm = ma* cj ; 

Ft = Fp + Fm; // Total thrust 

sfc=(ma/A_F)*3600/Ft; 

disp (” N (roundoff error)” ,F,” Total design thrust 

/s = ”); 

disp ( ”N (roundoff err or ) ” , Ft , ” Total thrust /s = ” 

) ; 

disp (”kg/ N thrust h”,sfc, ”Specific fuel 
consumption = ” ) ; 


Scilab code Exa 7.3 Calculation of the velocity leaving the nozzle 
clc ; 

p03=4.5; // Pressure at turbine inlet in bar 
T03=800 + 273; // Temperature at turbine inlet in 
k e 1 v i n 

p04=1.75; // Pressure at turbine outlet in bar 

eff_T = 0.75; //Turbine efficiency 
p05=1.03; // Pressure at state 5 in bar 
Cp = 1.05; // Specific lieat at constant pressure in 

kJ/kg K 

r = 1.38; // Specific lieat ratio 

T04=T03*(l-eff_T*(l-(l/(p03/p04)~((r-1)/r)))); // 

Temperature at state 4 

cj = sqrt (2*Cp*10~3*T04*(l-(l/(p04/p05)~((r-l)/r)))); 

// Yelocity 1 e a v i n g nozzle 
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disp (”K” , T04 , ” (i ) . Temperature of the gas entering 
the jet (nozzle) = ”); 

disp (”m/s” , cj , ” ( i i ) . Velocity of gas leaving the jet 
= ”); 


Scilab code Exa 7.4 Calculation of propulsive ethciency 
clc ; 

cj=2700; // The effective jet velocity from jet 

engine in m/s 

ci = 1350; // Flight velocity in m/s 

ma = 78.6; // Air flow rate in m/s 

a=ci/cj ; 

F = ma* ( c j - ci ) ; // Thrust 
P = F*ci; // Thrust power 

eff_P=2*a/(a + l) ; // Propulsive efficiency 

disp (”N” ,F,” ( i ) . Thrust = ” ) ; 
disp (”MN” ,P/10~6,” ( ii ) . Thrust power = ”); 
disp (”%” , ef f _P*100 , ” ( i i i ) . Propulsive efficiency = 
”); 


Scilab code Exa 7.5 Calculation of the thrust and Specihc fuel consump- 
tion 

clc ; 

pa = 0.458; // Ambient pressure in bar 

Ta = 248; // Ambient temperature in kelvin 

Ci=805*1000/3600; // Speed of the aircraft in m/s 

rp=4;// Pressure ratio 

DelP_comb =0.21 ; // Combustion chamber pressure loss 

in bar 
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T03 = 1100; // Turbine inlet temperature in k e 1 v i n 

eff_ram = 0.95; // Intake duct efficiency 

eff_c=0.85; // Compressor efficiency 

eff_T = 0.90; // Turbine efficiency 

eff_m = 0.99; // Mechanical efficiency of transmission 

eff_nozzle=0.95; // Nozzle efficiency 
CV = 43; // Low calorific value in MJ/kg 

Ac = 0.0935; // Nozzle outlet area in m' 2 

Cpa = 1.005;// Specific heat of air at constant 

pressure in kJ/kg K 

Cpg = 1.147;// Specific heat of fuel at constant 
pressure in kJ/kg K 
rg = 1.33;// Specific heat ratio of fuel 
r = 1.4; // Specific lieat ratio of air 

R = 287; // Characterist ic gas constant in J/kg K 

p01=pa*(l+eff_ram*((l+Ci~2/(2*Cpa*Ta*10~3))~(r/(r-1) 
)- 1 )) ; 
p02=p01*rp; 

T01 = Ta + Ci~2/(2*Cpa*10~3) ; 

T02=T01+T01*(rp~((r-1)/r)-1)/eff_c; 
T04=T03-(Cpa*(T02-T01))/(Cpg*eff_m); 
p03 = p02-DelP_comb ; 

T_04 = T03-(T03-T04)/eff_T ; 
p04 = p03*(T_04/T03)~(r/(r-l)) ; 

p04_pc=l/(l-(((rg-l)/(rg+1))/eff_nozzle))~(rg/(rg-1) 

) ; 

Tc=T04*(l/p04_pc)~((rg-1)/rg); 
pc=p04/p04_pc; 
row_c=(pc*10~5)/(R*Tc); 
cj=sqrt (rg*R*Tc); 
m = row_c*Ac*cj ; 

F=m*(cj-Ci)+Ac*(pc-pa)* 10~5; // Total thrust 

mf = (m*Cpg*(T03-T02))/(CV*10~3) ; 

sf c=mf *3600/F; // specific fuel consumption 

disp (”N (roundoff err or ) ” , F , ” Total thrust = ”); 
disp (”kg/N h (roundoff error)”,sfc,”specific fuel 
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consumption = ”); 


Scilab code Exa 7.6 Calculation of specific power output 
clc ; 

ci=600*1000/3600; // Velocity in m/s 
Cpa = 1.005;// Specific heat of air at constant 

pressure in kJ/kg K 

Cpg = 1.147;// Specific heat of fuel at constant 
pressure in kJ/kg K 
rg = 1.33;// Specific lieat ratio of fucl 
r = 1.4; // Specific lieat ratio of air 

R = 287; // Characteristic gas constant in J/kg K 

pa = 0.458; // Ambient pressure in bar 

Ta = -15 + 273; // Ambient temperature in kelvin 

rp = 9; // pressure ratio 

T03 = 1200; // Maximum temperature in ke 1 vin 

eff_ram = 0.9; // Intake duct efficiency 

eff_c=0.89; // Compressor efficiency 

eff_T = 0.93; // Turbine efficiency 

eff_m = 0.98; // Mechanical efficiency of transmission 

cj =ci 

T_01=Ta+(ci~2/(2*Cpa*10~3)) ; 
p_01=pa*(T_01/Ta)~(r/(r-l)); 
p01=eff_ram*(p_01-pa); 
p02=rp*p01; 

T01=T_01; 

T_02=T01*rp~((r-l)/r); 

T02=T01+(T_02-T01)/(eff_c); 

T_04=T03*(l/rp)~((rg-l)/rg); 

T04=T03-eff_T*(T03-T_04); 

WN=Cpg*(T03-T04)-Cpa*(T02-TO1)/eff_m; // net work 
done 

ef f _th = WN / ( Cpg* (T03-T02) ) ; // Thermal efficiency 
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disp (”kJ/kg (roundoff er r or ) ” , WN , ” Net work done = 
”); 

disp (”%” , ef f _th*100 , ” Thermal efficiency =”); 


Scilab code Exa 7.7 Determination of rate of fnel consumption 
clc ; 

pa = 0.7; // Ambient pressure in bar 

Ta = l + 273; // Ambient temperature in kelvin 

Ci=800*1000/3600 ; // Speed of the aircradt in m/s 
rp=5;// Pressure ratio 

eff_ram = 1.00; // Intake duct efficiency 

eff_c=0.85; // Compressor efficiency 

eff_T = 0.90; // Turbine efficiency 

ef f _comb=0.98 ; //Combustion efficiency 

eff_nozzle=0.95; // Nozzle efficiency 
rp_T=2.23;// Turbine pressure ratio 
CV = 43; // Low calorific value in MJ/kg 
Cpa = 1.005;// Specific heat of air at constant 
pressure in kJ/kg K 

Cpg = 1.005;// Specific heat of fuel at constant 
pressure in kJ/kg K 
rg = 1.4;// Specific lieat ratio of fuel 
r = 1.4; // Specific lieat ratio of air 

R = 287; // Characteristic gas constant in J/kg K 

F = 25000 ; // Thrust in N 

c j=2*Ci ; 

T_01=Ta+(Ci~2/(2*Cpa*10~3)) ; 

T01=T_01; 

T02=T01+(T01*(((rp)~((r-l)/r))-l))/eff_c; 
p_01=pa*(l+Ci~2/(2*Cpa*10~3*Ta))~(r/(r-l)); 
p01=eff_ram*(p_01-pa); 
p02=rp*p01; 
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T03 = (T02-T01)/(ef f _T *(1-1/rp_T~((r-1)/r))); 
ma=F/(cj-Ci); 

// Neglecting the effect of the mass addition of 
fuel on the riglit hand side 

mf = (ma*Cpa*(T03-T02))/(eff_comb*C¥*10"3) ; 

disp (” kg/s ” , ma , ” Mass flow rate of air = ”); 
disp (”kg/s (roundoff er ror ) ” , mf , ” Mass flow rate 
o f fuel = ”) ; 


Scilab code Exa 7.8 Calculation of the take off thrust 
clc ; 

Ta = 288; // Ambient temperature in kelvin 

pa = 1.01; // Ambient pressure in bar 

p04=2.4; // Stagnation pressure in bar 

T04 = 1000;// Stagnation temperature in kelvin 
m = 23; // Mass flow rate in kg/s 

rp=1.75; // Pressure ratio 
eff_f=0.88 ; // Efficiency of the fan 
eff_ft=0.9; // Efficiency of the fan turbine 

Cpa = 1.005;// Specific heat of air at constant 

pressure in kJ/kg K 

Cpg = 1.147;// Specific lieat of fuel at constant 
pressure in kJ/kg K 
rg = 1.33;// Specific heat ratio of fuel 
r = 1.4; // Specific lieat ratio of air 

R = 284.6 ; // Characteristic gas constant in J/kg K 

T01=Ta; 
pO1=pa; 

pc = p04*(2/(r + 1))~(r/(r-l)) ; 

// since pc>pa the nozzle will choke 
Tc=T04*(2/(r+1)); 
row_c = pc*10~5/(R*Tc) ; 
cj = sqrt (r*R*Tc) ; 


53 



22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

1 

2 

3 

4 

5 

6 

7 

8 


A=m/(row_c * cj) ; 
pl=pa; 

F=m*cj + (A*(pc-pl)* 10~5) ; 

// For fan engine 

T_02 = T01*(rp)~((r-l)/r) ; 

T02 = T01 + (T_02-T01)/eff_f ; 

// For cold nozzle 

m_nozzle=2*m; // Flow through cold nozzlc 

pcl=p01*rp*(2/(r+1))“(r/(r-1)); 

F_cold = m_nozzle* sqrt (2*Cpa* 10~3*(T02-T01)) ; 

// Fan Turbine 

T05 = T04-((m_nozzle*Cpa*(T02-T01))/(m*Cpg)) ; 
T_05=T04-(T04-T05)/eff_ft; 
p_05 = p04*(T_05/T04)~(rg/(rg-l) ) ; 
pc=p_05*(2/(rg+1))~(rg/(rg-l)); 

F_hot=m*sqrt (2*Cpg* 10~3*(T05-TO1)) ; 

Takeoffthrust= F_cold + F_hot; 

disp (”m~2 (roundoff error ) ” , A , ” Nozzle Exit area = 

disp (”N (roundoff error ) ” , F , ” Total Thrust = ”); 

disp (”N (roundoff error ) ” , Takeof f thrust , ” Take-off 

Tlirust = ” ) ; 


Scilab code Exa 7.9 Calculation of thrust provided by the engine and the 
thrust power developed 

clc ; 

ma = 18.2; // Massflow rater in m/s 

Mi = 0.6 ; // Mach number 

pa = 0.55; // Ambient pressure in bar 

Ta = 255; // Ambient temperature in kelvin 

rp = 5; // Pressure ratio 

T03 = 1273; // Maximum temperature in kelvin 

eff_c=0.81; // Compressor efficiency 
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eff_T = 0.85; // Turbine efficiency 

eff_nozzle=0.915; // Nozzle efficiency 
ef f _ram = 0.9 ; // Intake dnct efficiency 

CV=45870; // Low calorific value in kJ/kg 
Cpa = 1.005;// Specific lieat of air at constant 
pressure in kJ/kg K 

Cpg = 1.147;// Specific heat of fuel at constant 
pressure in kJ/kg K 
rg = 1.33;// Specific lieat ratio of fuel 
r = 1.4; // Specific heat ratio of air 

R = 284.6; // Characteristic gas constant in J/kg K 

ci=Mi* sqrt (r*R*Ta); 

T_01=Ta+ci~2/(2*Cpa*10~3); 

T01=T_01; 

p_01=pa*(T01/Ta)~(r/(r-01)) ; 
p01=eff_ram*(p_01-pa)+pa; 
p02=rp*p01; 

T02=T01*(l+((rp~((r-l)/r))-l)/eff_c); 

Wc=ma*Cpa*(T02-T01); 

WT=Wc; 

mf=ma/((CV/(Cpg*(T03-T02)))-1); 
f1=mf/ma; 

T04=T03-(WT/((ma+mf)*Cpg)); 

rp_T=(l/(l-((l-(T04/T03))/eff_T)))~(r/(r-l)); 
p03=p02; 
p04=p03/rp_T; 

p04_pc=l/(l-((rg-l)/((rg+l)*eff_nozzle)))~(rg/(rg-l) 

) ; 

pc=p04_pc/p04; 

Tc=T04*(l/p04_pc)~((rg-l)/rg); 
cj = sqrt (r*R*Tc) ; 
row_c=pc*10~5/(R*Tc); 

An=(ma+mf)/(row_c*cj); 

F=(ma+mf)*cj-ma*ci+An*(pc-pa); 

Fp = F* ci ; 

disp (”kW ( roundoff error ) ” , Wc , ” Work of 
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compression = ” ) ; 

disp (”kW (roundoff err or ) ” , WT , ” Power output of 
t h e t u r b i n e = ” ) ; 
disp (f 1 , ” Fuel— Air ratio = ”); 
disp (”N (roundoff error)” , F , ” Thrust = ”); 
disp (”kW ( roundoff error ) ” , Fp/1000 , ” Thrust power 
= ”); 


Scilab code Exa 7.10 Calculation of exit speed of the gases and the thrust 
developed 

clc ; 

ma= ( 12* 10 ~ 4)/3600 ; // Air flow rate in kg/s 

T01 = 15 + 273; // Temperature in kelvin 

rp=4; // pressure ratio 

p01=1.03; // Pressure in bar 
T02 = 182 + 273; // Temperature in kelvin 
T03=815 + 273; // Temperature in kelvin 

T04 = 650 + 273; // Temperature in kelvin 

ci=800*1000/3600; // Velocity in m/s 
eff_nozzle=0.90; // Nozzle efficiency 
Cpa = 1.005;// Specific heat of air at constant 
pressure in kJ/kg K 

Cpg = 1.147;// Specific lieat of fuel at constant 
pressure in kJ/kg K 
rg = 1.33;// Specific heat ratio of fuel 
r = 1.4; // Specific lieat ratio of air 

p03 =4.12; // in bar 

eff_c=l/((T02-T01)/(T01*((rp~((r-l)/r))-l))); 
eff_T = eff_c ; 

Wc = ma*Cpa*(T02-T01) ; 

rp_T = (l/(l-((T03-T04)/(eff_T*T03))))~((r/(r-l))) ; 
p04=p03/rp_T; 

p04_pc = l/(1 -((rg-1)/((rg + 1)*eff_nozzle)))~(rg/(rg-1) 
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) ; 

p5=p01; 

T_5=T04*(p5/p04)~((rg-l)/rg); 
T5 = T04-eff_nozzle*(T04-T_5) ; 
cj =sqrt (2*Cpg*10~3*(T04-T5)); 
F=ma*(cj-ci); 


disp (”%” , ef f _c*100 , ” Efficiency of the 

”) ; 

disp (”%” , ef f _T*100 , ” Efficiency of the 
disp (”kW” , Wc , ” Compressor work = ”); 
disp (”m/s (roundoff error ) ” , cj , ” The 
gases = ”); 

disp (”N (roundoff err or ) ” , F , ” Thrust 
) ; 


compressor = 
Turbine = ” ) ; 
exit speed of 
developed = 
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Chapter 8 

Centrifugal Compressors 


Scilab code Exa 8.1 Calculation of compressor efficiency 
clc ; 

N = 11500; // Speed in rpm 

T01=21 + 273; // Inlet total temperature in kel vin 

p01 = 1; // Inlet total pressiire in bar 
p02=4;// Outlet total pressure in bar 
D = 0.7 5 ; / / i m p e 11 e r d i a m e t e r i n m 
mu=0.92;// slip factor 

Cp = 1.005; // specific heat at constant pressure in 
kJ/kg K 

r = 1.4; // Specific lieat ratio 

u = 3.14 + D + N/60 ; 

W=mu*u~2; 

T02=W/(Cp*10~3)+T01; 

T_02=T01*(p02/p01)“((r-1)/r); 
eff_c = (T_02-T01)/(T02-T01) ; 

disp (”%”,eff_c*100,”Efficiency of the compressor = 

”); 
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Scilab code Exa 8.2 Estimation of the probablc axial widtli of the impellcr 


clc ; 

m = 35; // mass flow rate of air in kg/s 

D = 0.7 6 ; / / I m p e 11 e r d i a m e t e r i n m 
N = 11500; // speed in rpm 

eff_c=0.8; // Efficiency of the compressor 

rp=4.2; // Pressure ratio 

cr = 120; // Radial velocity in m/s 

p 01 = 1; // Inlet pressure in bar 

T01=47 + 273; // Inlet temperature in kelvin 

Cp = 1.005; // specific lieat at constant pressure in 

kJ/kg K 

r = 1.4; // Specific lieat ratio 

R = 287; // Characteristic gas constant in J/kg K 

T_02 = T01*rp~((r-l)/r) ; 

T02 = T01 + (T_02-T01)/eff_c ; 

// ignoring the effects of the velocity of flow 

p02=rp/pO1; 

row2=p02*10~5/(R*T02); 

Atip=m/(row2*cr) ; 

AW = Atip/(3 . 14*D) ; // Axial width 
disp (”cm” , AW*100 , ” Axial Widtli = ”); 


Scilab code Exa 8.3 Calculation of theoretical blade angles 
clc ; 

D = 0.15; // Inlet eye diameter in m 

N = 20000; // Speed in rpm 

cal = 107; // Axial velocity in m/s 
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T01=294; // Inlet temperature in k e 1 v i n 

p01 = 1.03; // Inlet pressure in kg/crn' 2 

Cp = 1.005; // specific lieat at constant pressure in 
kJ/kg K 

r = 1.4; // Specific lieat ratio 

R = 287; // Characteristic gas constant in J/kg K 

ul=3.14*D*N/60; 

beta_l=atand (cal/ul);// Blade angle 

cr=ul/cosd (beta_l); 

a=sqrt (r*R*(T01-cal~2/(2*Cp*10~3))); 

M=cr/a; // Mach number at the tip 

disp (” degree” , beta_l , ” ( i ) . Theoretical angle of the 
blade at this point = ”); 

disp (M, ” ( i i ) . Mach number of the flow at the tip of 
the eye = ”); 


Scilab code Exa 8.4 Calculation of hnal temperature of the gases and the 
work done per kg of gas 

clc ; 

T01=0 + 273; // Inlet gas temperature in k e 1 v i n 

p01=0.7; // Inlet pressure in bar 

p02 = 1.05; // Delivery pressiire in bar 
eff_c=0.83; // Compressor efficiency 

Cp = 1.005;// Specific lieat at constant pressure in 

kJ/kg K 

Cv = 0.717;// Specific heat at constant volume in kJ 

/kg K 

r = 1.4; // Specific heat ratio 

T_02 = T01*(p02/p01)“((r-l)/r) ; 

T02 = T01 + (T_02-TOl)/ef f _c ; // Final temperature of 
the gas 
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Wc = Cp* (T02-T01) ; // Work of compression 

// With additional compressor 

T_03 = T02*(p02/p01)~((r-l)/r) ; 

T03 = T02 + (T_03-T02)/eff_c ; 

T_03 = T01*(p02/p01)~(2*(r-1)/r) ; 
eff_overall=(T_03-T01)/(T03-T01); 

disp (”K” , T02 , ” Final temperature of the gas = ”); 

disp (” k J/kg” , Wc , ” Work of compression = ”); 

disp (”%” , eff_overall*100 ,” Overall efficiency =”); 


Scilab code Exa 8.5 Calculation of impeller diameters and the width at 
the impcller exit 

clc ; 

N = 12500; // Speed in rpm 

ra = 15; // Mass flow rate in kg/s 

rp=4; // Pressure ratio 

eff_c=0.75; // Isentropic efficiency 

mu = 0.9; // Slip factor 

pi=0.3; // Flow coefficient at impeller exit 

D = 0.15; // Hub diameter in m 

ca2 = 150; // Axial velocity in m/s 

T01=275; // Inlet temperature in kelvin 

p 01 = 1; // Inlet pressure in bar 

Cp = 1.005;// Specific lieat at constant pressure in 

kJ/kg K 

Cv = 0.717;// Specific heat at constant volume in kJ 

/kg K 

r = 1.4; // Specific heat ratio 

R = 287; // Characteristic gas constant in J/kg K 

u2=ca2/pi; 

P=m*mu*u2~2/1000; // Power output 
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D2 = u2 *60/(3.14*N) ; 

Tl = T01-ca2~2/(2*Cp*10~3) ; 
pl=p01*(Tl/T01)~(r/(r-1)); 
rowl=pl*10~5/(R*T1); 

Al=m/(rowl*ca2); 

D1 = sqrt ((A1*4/(3.14))+D~2) ; 
p3_pl=rp; 

p2=2*pl; 

T_2=Tl*(p2/pl)~((r-l)/r); 

T2=T1+(T_2-T1)/eff_c; 
row2=p2*10~5/(R*T2); 

W2=(m)/(row2*ca2*3.14*D2); 

disp (”kW” ,P ,”Power = ” ) ; 
disp (” Impe 11 er Diameters”); 

disp (”cm” , D2*100 , ”D2 = ” , ” cm (roundoff error)”,Dl 
*100 ,”D1 = ” ) ; 
d i s p (” I m p e 11 e r w i d t h ” ) 

disp (”cm ( roundoff error ) ” , W2*100 , ”W2 = ” ) ; 


Scilab code Exa 8.6 Calculation of the minimum possible depth of the 
diffuser 

clc ; 

m = 14; // mass flow rate in kg/s 

rp=4; // pressure ratio 

N = 12000; // Speed in rpm 

T01=288; // Inlet temperature in k e 1 v i n 

p01=1.033; // Inlet pressure in bar 

Cp = 1.005;// Specific lieat at constant pressure in 

kJ/kg K 

Cv = 0.717;// Specific heat at constant volume in kJ 

/kg K 

r = 1.4; // Specific lieat ratio 

R = 287; // Characterist ic gas constant in J/kg K 
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mu = 0.9; // Slip factor 

chi=1.04; // Power input factor 

eff_c=0.8; // Compressor efficiency 

T03=(((rp~((r-1)/r))-1)*T01/eff_c)+T01;; 

U=sqrt ((T03-T01)*Cp* 10~3/(chi*mu)) ; 

D=U *60/(3.14*N) ; 

T3=T03/1.2; 
c2=sqrt (r*R*T3); 
ca2=sqrt (c2~2-(mu*U)~2); 

T02=eff_c*(T03-T01)+T01; 

Loss =T03-T02; 

T2=T3-Loss/2 

p2 = p01*(T2/T01)~(r/(r-1)) ; 
row2=p2*10~5/(R*T2); 

A=m/(row2 * ca2) ; 

Depth = A/(2*3.14 + D/2) ; 

disp (”cm” , D*100 , ” O verall diameter of the Impeller 
= ”) ; 

disp (”cm ( roundoff error ) ” , Depth*100 , ” Depth of 
t h e d i f f u s e r = ” ) ; 


Scilab code Exa 8.7 Calculation of impeller and diffuser blade angles at 
inlet 

clc ; 

N = 10000; // Speed in rpm 

Q = 600; // Flow rate m~2/min 

rp=4; // Pressure ratio 

eff_c=0.82; // Compressor efficiency 

T01=293; // Inlet temperature in k e 1 v i n 

pO1=1.0; // Inlet pressure in bar 

Cp = 1.005;// Specific lieat at constant pressure in 
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kJ/kg K 

Cv = 0.717;// Specific lieat at constant volnme in kJ 
/kg K 

r = 1.4; // Specific lieat ratio 

R = 287; // Characteristic gas constant in J/kg K 

ca = 60; // Axial velocity in m/s 

D2_D1=2 ;// Diameter ratio 

T_03 = T01*rp"((r-1)/r) ; 

T03 = T01 + (T_03-T01)/eff_c ; 
u2 = sqrt (Cp* 10~3*(T03-TO1)) ; 

Wc = u2~2; // Work of compression 

D2 = (u2*60/(3.14*N)) ; 

D1=D2/D2_D1; 

T1 = T01-(ca~2/(2-Cp*10~3)) ; 
pl=p01*(Tl/T01)~(r/(r-1)); 
rowl=pl*10~5/(R*T1); 

Wroot = (Q/60)*(1/(ca*3.14* D1)) ; 
ul=3.14*N*D1/60; 
alpha_root=atand (ca/ul); 
alpha_tip= atand (ca/u2); 

disp (”(i).Power input ”); 
disp (”kW/kg/ s” ,Wc/1000 , ”Wc = ”); 
disp (” ( i i ) . Impe 11 er Diameters”); 
disp (”m”,D2,”D2 = ”,”m”,Dl,”Dl = ”); 
disp (” ( i i i ) . Impeller and diffuser blade angles at 
i n 1 e t ” ) ; 

disp (” degree” , alpha_tip , ” alpha_tip = ”,”degree”, 
alpha_root,”alpha_root = ”); 


Scilab code Exa 8.8 Calculation of slip factor 
clc ; 

rp=4; // Pressure ratio 
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eff_c=0.8; // Compressor efficiency 

N = 15000; // Speed in rpm 

T01=293; // Inlet temperatnre in kelvin 

De=0.25; // Diameter of eye in m 

Cl = 150; // Absolute velocity in m/s 

D i = 0.6 ; / / I m p e 11 e r d i a m e t e r i n m 

al=25; // in degree 

Cp = 1.005;// Specific heat at constant pressure in 

kJ/kg K 

Cv = 0.717;// Specific heat at constant volume in kJ 

/kg K 

r = 1.4; // Specific lieat ratio 

R = 287; // Characterist ic gas constant in J/kg K 

T02 = T01*rp~((r-l)/r) ; 

DelT_actual = (T02-T01)/eff_c ; 

P=Cp*10~3*DelT_actual ; // Power input 
ul=3.14*De*N/60; 
ctl=Cl*sind (al); 

// At Exit 
u2 =3.14*Di*N/60; 
ct2=(P+(ul*ctl))/u2; 
mu=ct2/u2; // Slip factor 

disp (mu Slip Factor = ”); 


Scilab code Exa 8.9 Determination of number of radial impeller vanes 
clc ; 

P = 180*10~3; // Power input in J 
N = 15000; // Speed in rpm 
al=25; // in degrees 

De=0.25; // Mean dia of the eye in m 

Di = 0.6 ; // Impeller tip diameter in m 
cl = 150; // Absolute air velocity at inlet in m/s 
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ul=3.14*De*N/60 ; 
u2 =3.14*Di*N/60; 
ctl=cl*sind (al); 
ct2=(P+(ul*ctl))/u2; 
mu= ct2/u2; 

z=(1.98)/(l-mu); // Number of impeller vanes 
di sp (z , ” Nurnber of impeller vanes using Stanitz 
formulae = ” ) ; 


Scilab code Exa 8.10 Calculation of the torque power required and the 
head developed 

clc ; 

m = 30; // mass flow rate in kg/s 

N = 15000; // Speed in rpm 

r2=0.3; // Radius in m 

D2 = r2*2; // Diameter in m 

w2 = 100; // Relative velocity in m/s 

beta_l=80; // in degrees 
p 01 = 1; // Inlet pressure in bar 

T01=300 // Inlet temperature in k e 1 v i n 
Cp = 1.005; // specific lieat at constant pressure in 

kJ/kg K 

r = 1.4; // Specific lieat ratio 

R = 287; // Characteristic gas constant in J/kg K 

u2 =3.14*D2*N/60; 

ct2=u2-(w2*cosd (beta_l)); 

Fr=m*ct2*r2; 

P = Fr*(2*3.14*N/60) ; 

W = u2 * ct2; 

P02=p01*(l+(W*10~-3/(Cp*T01)))~(r/(r-l)); 
disp (”Nm” ,Fr,”Torque = ”); 
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22 disp (”kW” , P/1000 , ” Power = ”); 

23 disp (” bar ” , P02 , ” Head Developed = ”); 
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Chapter 9 

Axial Flow Compressors 


Scilab code Exa 9.1 Estimation of blade angles 
clc ; 

n = 10; // No of stages in the axial flow compressor 

rp = 5; // Overall pressure ratio 

eff_C = 0.87; // Overall isentropic efficiency 

Tl = 15 + 273; // Temperature of air at inlet in kelvin 

u = 210; // Blade speed in m/s 

ca = 170; // Axial velocity in m/s 

wf = l; // Work factor 

r = 1.33; // Specific lieat ratio 

Cp = 1.005; // Specific lieat in kJ/kg K 

Del_Tstage=(Tl*(rp~((r-l)/r)-1))/(n*eff_C); // 

Temperature increase per stage 
// By property relations and let us assume 
// tan_betal—tan_beta2 = Del_Tstage*Cp/(wf*u*ca) 

// tan_betal + tan_beta2=u/ca for 50% reaction 
// To solve this above equations using matrix rnethod 
a=[l,-1;1 ,1] ; c=[(Del_Tstage*Cp*10~3/(wf*u*ca)) ;u/ca 
] ; 

b=a\c; 

betal=atand(b(1)); // Blade angles at inlet 
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beta2 = atand (b (2) ) ; // Blade angles at outlct 


disp (” degree (roundoff error ) ” , beta2 , ” Blade angle 
at outlet = ”,”degree ( roundoff e r r or ) ” , betal , 
” Blade angle at inlet = ”); 


Scilab code Exa 9.2 Calculation of rnass flow rate and degree of reaction 
clc ; 

Pl=1.0132; // Inlct air pressure in bar 

T01=288; // Inlet air temperature in k e 1 v i n 

ca=150; // axial velocity in m/s 

dtip = 60; // Tip diameter of rotor in crn 
dhub=50; // Hub diameter of rotor in crn 
N=100; // Speed of rotor in rps 

t_angle=30; // Deflected angle of air in degree ( in 
question it is 30.2 but in solution it is 30) 
P2_P1=1.2; // Stage pressure ratio 

Cp = 1005; // Specific heat in J/kg K 
r = 1.4; // Specific heat ratio 

R = 287; // Characteristic gas constant in J/kg K 

u=(3.142857*(dhub+dtip)*10~-2*N)/2; // Mean blade 

velocity 

beta_l=atand(u/ca); // Blade angle at inlet 
beta_2=beta_l-t_angle; // As air is deflected by 30 
// from velocity triangle 
x=ca*tand(beta_2); 
alpha_2=atand ((u-x)/ca); 

Cl=ca; 

T1 = T01-(Cl “2/(2*Cp) ) ; // Static temperature at inlet 
P2=P1*P2_P1; // Pressure at outlet 

T2 = T1 * (P2/P1) “ ((r-1)/r) ; // Static temperature at 
outlet 

row_2 = (P2 * 10~5) / (R*T2) ; // Density at outlet 
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ra = 3.14* ( dt ip "2-dhub ~2) * ca* row_2 * 10~-4/4 ; // Mass 

flow rate 

wf = l; // Work factor 

P=wf*u*ca*m*(tand(beta_l)-tand(beta_2))/1000; // 

Power developed 

R=ca*(tand(beta_1)+tand(beta_2))/(2*u); // Degree of 
r e a c t i o n 

disp (” kg/s ” , m , ” Mass flow rate = ”); 
disp(”kW ( Error due to more decimal values in 
expression ) ” , P , ” Power devcloped = ”); 
disp (R,”Degree of Reaction = ”); 


Scilab code Exa 9.3 Estimation of number of stages of the compressors 
clc ; 

beta_l=45; // Inlet blade angle in degree 

beta_2 = 10; // Outlet blade angle in degree 

rp = 6; // Compressor pressure ratio 

eff_C = 0.85; // Overall isentropic efficiency 

Tl=37 + 273; // Inet static temperature in kelvin 

u = 200; // Blade speed in m/s 

Cp = 1005; // Specific lieat in J/kg K 

r = 1.4; // Specific lieat ratio 

R = 287; // Characteristic gas constant in J/kg K 

// ( i ) . wf=l 

wf = l; // Work factor 

ca = u/ ( tand ( beta_ 1) +tand (beta_2 ) ) ; // Axial velocity 

Del_Tstage=wf*u*ca*(tand(beta_l)-tand(beta_2))/Cp; 

// Stage temperature drop 

Del_Toverall=(Tl*(rp~((r-l)/r)-l))/eff_C; // Overall 
temperature drop 

n=Del_Toverall/Del_Tstage ; // No of stages 
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disp (n,”Number of stages required = ” , ” (i ) . wf = 1”) 


// (ii).wf = 0.87 
wf =0.87; // Work factor 

ca = u/(tand (beta_ 1)+tand (beta_2 ) ) ; // Axial velocity 

Del_Tstage=wf*u*ca*(tand(beta_l)-tand(beta_2))/Cp; 

// Stage temperature drop 

Del_Toverall=Tl*(rp"((r-l)/r)-l)/eff_C; // Overall 
temperature drop 

n=Del_Toverall/Del_Tstage; // No of stages 

disp (n,”Number of stages required = ” , ” ( i i ) . wf = 
0.87”); 


Scilab code Exa 9.4 Determination of Mach number relative to Rotor 
clc ; 

rp=4; // Total head pressure ratio 

eff_0=0.85; // Overall total head isentropic 

efficiency 

T01=290; // Total head inlet temperature in kelvin 

alpha_l=10; // Inlet air angle in degree 
alpha_2=45; // Outlet air angle in degree 

u = 220; // Blade velocity in m/s 

wf=0.86; // Wok done factor 

R = 284.6; // Characteristic gas constant in kJ/kg K 

Cp = 1005; // Specific heat in J/kg K 
r = 1.4; // Specific heat ratio 

eff_P=l/(logl0(((rp~((r-l)/r)-l)/eff_0)+l)/(logl0(rp 
)*((r-1)/r)));; 

// From velocity triangle 

ca=u/ (tand(alpha_1)+tand(alpha_2)); // Axial 

velocity 
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Del_Tstage=wf*u*ca*(tand(alpha_2)-tand(alpha_l))/Cp; 

// Stage temperature rise 

T02 = T01*(rp)“((r-l)/(r*eff_P)) ; // Total head 
temperature 

T02_T01 =T02-T01 ; // Total temperature rise 
n = T02_T01/Del_Tstage ; // Total number of stages 
// from velocty traingles 

wl=ca/cosd(alpha_2); 
c1=ca/cosd(alpha_1); 

T1 = T01 - c 1 ~ 2/(2* Cp) ; // Static temperature 
M = wl /sqrt (r*R*Tl) ; // Mach number at inlet 

disp (eff_P*100,”Polytropic efficiency of the 
compressor = ”); 

disp (n,”Total number of stages = ”); 
disp (M,”Mach number at inlet = ”); 


Scilab code Exa 9.5 Calculation of pressure rise per blade ring and the 
power input per stage 

clc ; 

Q = 1000; // Flow rate of free air in nP3/min 

Pl=0.98; // Inlet pressure in bar 

Tl = 15 + 273; // Inlet temperature in kelvin 
Dm = 0.6; // Mean diameter in m 

h = 6.75; // blade lengtli in cm 

CL=0.6; CD=0.05; // At zero angle of incidence 

Cp = 1.005; // Specific heat in kJ/kg K 
r = 1.4; // Specific heat ratio 

R = 287; // Characteristic gas constant in J/kg K 
k = l-0.1; //Blade occupys 10% of axial area 
N = 6000; // speed in rpm 

Ac = 19.25*10~-4; // Projected area in m' 2 

n = 50 ; 

eff_C = l; // Efficiency of compressor 
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row = (Pl*10~5)/(R*Tl); // Density 

A=k*3.14*Dm*h*10~-2; // Area of axial 

ca = Q/(60*A); // Axial velocity 

u=3.14*Dm*N/60; // Blade velocity 

beta_l=atand(u/ca); // Blade angle at inlet 

w=sqrt (ca~2 + u~2); // From velocity triangle 

L=CL*row*w~2*Ac/2; 

D=CD*row*w~2*Ac/2; 

P=(L*cosd(beta_1)+D*sind (beta_l))*u*n*10~-3; // 

Power input / stage 
m = Q*row/60 ; // mass flow rate 

rp=((P*eff_C/(m*Cp*Tl))+l)~(r/(r-l)); // pressure 
r a t i o 

P2 = rp*Pl ; // Pressure 

disp (”kW (Roundoff error )”,P,” Power input/stage 
= ”) ; 

disp (” bar” , P2 , ” Pressure at outlet = ”); 


Scilab code Exa 9.6 Determination of the direction of the air entry to and 
exit from the rotor 

clc ; 

Tl=290; // Temperature at inlet in kelvin 

n = 10; // Number of stages 

rp=6.5; // Pressure ratio 

m = 3; // mass flow rate in kg/s 

eff_C=0.9; // isentropic efficiency of the 

compression 

ca = 110; // Axial velocity in m/s 

u = 180; // Mean blade velocity in m/s 

Cp = 1.005; // Specific lieat in kJ/kg K 
r = 1.4; // Specific lieat ratio 

R = 287; // Characterist ic gas constant in J/kg K 
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13 T_2 = (rp) ~ ( (r-1)/r) *T1 ; // temperature after 

isentropic compression 

14 T2 = ((T_2-T1) /ef f _C)+T1 ; // Temperatnre after actnal 

compression 

15 P=m*Cp* (T2-T1) ; // Power given to the air 

16 Del_Tst age = (T2-T1)/n ; // Temperature rise per stage 

17 Del_ct = Cp* 10~3*Del_Tstage/u; // For work done per kg 

of air per second 

18 // To find blade angles let solve the following 

equations 

19 // Del_ct=ca (tan beta_l—tan beta_2) for symmetrical 

st ages 

20 // u=ca(tan beta_l=tan beta_2) for degree of 

reaction = 0.5 

21 // Solving by matrix method 

22 A=[l,—1;1,1] ; C= [Del_ct/ca;u/ca] ; 

23 B = A\C; 

24 // Blade angles at entry and exit 

25 beta_ 1 = atand (B (1) ) ; 

26 beta_2 = atand (B (2) ) ; 

27 

28 disp (”kW (roundoff er r or ) ” , P , ” Power given to the 

air = ”) ; 

29 disp (” degree” , beta_2 , ” Blade angle at exit = ”,” 

degree”,beta_1 , ” Blade angle at inlet = ”); 


Scilab code Exa 9.7 Calculation of the rotational speed and thc lcngth of 
the last stage 

1 clc ; 

2 rp=4; // Overall pressure ratio 

3 m = 3; // mass flow rate in kg/s 

4 eff_pc=0.88; // Polytropic efficiency 

5 Del_Tstage =25; // The stagnation temperature 
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pressure rise in kelvin 

6 cl = 165; // Absolnte velocity in m/s 

7 alpha_l=20; // air angle from axial direction in 

degree 

8 wf=0.83; // Workdone factor 

9 D = 18; // Mean diameter of the last stage rotor in cm 

10 P01 = 1.01; // Ambient pressure in bar 

11 T01=288; // Ambient temperature in kclvin 

12 Cp = 1005; // Specific lieat in J/kg K 

13 r = 1.4; // Specific heat ratio 

14 R = 287; // Characteristic gas constant in J/kg K 

15 

16 n = l/(1 - (r-1) / (r* ef f _pc ) ) ; 

17 T02 = T01 * (rp) “ ((n-1)/n) ; // Total pressure at stage 2 

18 Del_Toverall= T02-T01; // Overall temperature 

difference 

19 Ns = Del_Toverall/Del_Tstage ; // Number of stages 

20 eff_C=((rp“((r-l)/r)-l)/(rp~((r-l)/(r*eff_pc))-l)); 

// Efficiency of compressor 

21 rpl=(1+(eff_C*Del_Tstage/T01))“(r/(r-1)); // 

Pressure ratio acrocc first stage 

22 Del_Tstage 1 = Del_Toveral 1 /Ns ; // Temperature rise 

across stage 1 

23 T01 s =T02 -Del_Tstage 1 ; // Temperature at inlet to 

last stage 

24 rpls=(1+(eff_C*Del_Tstage1/TOls))~(r/(r-1)); // 

Pressure ratio acrocc last stage 

25 // For symmetrical blade , R=0.5 

26 beta_2 = alpha_l ; 

27 ca=cl*cosd (alpha_l); // Axial velocity 

28 beta_l=atand(sqrt (((Cp*Del_Tstagel/(wf*ca))/ca)+( 

tand(beta_2))~2)); // blade angle 

29 u= ca* (t and (bet a_ 1) +t and (bet a_2 ) ) ; // mean velocity 

of blade 

30 N = 60*u/(3.14*D*10~-2*60) ; // Speed in rps 

31 Po=rp/rpls; // Total pressure at inlet to the last 

st age 

32 T0 = T01s ; // Total temperature to the last stage 
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Tst =T0 - c 1 ~ 2/(2* Cp) ; // Static temperature 
Pst=Po/(TO/Tst)“((r-1)/r); // Static pressure 
row = (Pst*10~5)/(R*Tst); // Density 

h=m/(ca*row*3.14*D*10~-2); // Length of last stage 
disp (Ns Number of stages = ”); 

disp (rpl , ” Pressure ratio across first stage = ”); 
disp (” (roundoff error ) ” , rpls , ” Temperature at 
inlet to last stage = ”); 

disp (” degree (roundoff error ) ” , beta_l , ” betal=” ); 
disp (”rps (roundoff er ror ) ” , N , ” Speed = ”); 
disp (”cm (roundoff err or ) ” , h* 100 , ” Length of last 
stage = ”); 


Scilab code Exa 9.8 Calculation of the stage stagnation pressure ratio and 
the power input 

clc ; 

N = 6000; // Speed in rpm 

Del_rise=20; // Stagnation temperature rise in 
k e 1 v i n 

wf=0.93; // Work done factor eff_c=0.89; // 

Isentropic efficiency of the state 
cl = 140; // Inlet velocity in m/s 

p01 = 1 . 01 ; // Ambient pressure in bar 

T01=288; // Ambient temperature in kelvin 

Ml=0.95; // Mach number 

Cp = 1.005; // Specific heat in J/kg K 
r = 1.4; // Specific heat ratio 

R = 287; // Characteristic gas constant in J/kg K 

H_T_ratio=0.6; // Hub tip ratio in 
eff_s=0.89; // Stage efficiency 

Tl = T01-cl~2/(2*Cp*10~3) ; 
wl=Ml*sqrt (r*R*Tl); 
beta_l=acosd (cl/wl); 
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u=wl*sind (beta_l); 

beta_2=atand (tand(beta_l)-((Cp*10~3*Del_rise)/(u*wf 
*cl))); 

pl=p01/(TOl/Tl)~(r/(r-1)); 
row_l=(pl*10~5)/(R*T1); 

Rtip=60*u/(2*3.14*N); 

Rroot=H_T_ratio*Rtip; 

Rm=(Rtip+Rroot)/2; 
h=Rtip-Rroot; 
m=row_l*2*3.14*Rm*h*cl; 

rp=(l+(eff_s*Del_rise)/(T01))~(r/(r-l)); 
P=m*Cp*Del_rise; 
uroot=2*3.14*Rroot*N/60; 
beta_lroot=atand (uroot/cl); 

beta_2root = atand (tand (beta_lroot)-((Cp* 10~3* 
Del_rise)/(wf*uroot*cl))); 

disp (” degree” ,beta_2,”beta 2 = ” , ” degree” ,beta_l , ” 
beta 1 = ”,”Rotor air angles at tip : ” ,”m” , Rtip , ” 
Tip Radius = ” ,” ( i ) . ”) ; 

disp (”kg/s (Roundoff error ) ” ,m, ” Mass flow rate = 

” ( i i ) - ”) ; 

disp (”kW” , P , ” Power input = ” , rp , ” Stagnation 
pressure ratio = ” , ” ( i i i ) . ” ) ; 
disp (” degree” , beta_2root , ” beta 2 = ”,”degree”, 

beta_ lroot , ” beta 1 = ,: ,”Rotor air angles at root 
s e c t i o n s ” , ” ( i v ) . ” ) ; 


Scilab code Exa 9.9 Determination of the stage efficiency and the work 
done factor of an axial ffow compressor 

clc ; 

rp=1.35; // Actual pressure ratio 

DelT_rise=30; // Actual temperature rise in K 
beta_l=47; // Inlet blade angle in degree 
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beta_2 = 15; // Outlet blade angle in degree 

u = 225; // Peripheral velocity in m/s 

ca = 180; // Axial velocity in m/s 

T01 =27 + 273; // Ambient temperature in kelvin 

Cp = 1.005; // Specitic lieat in KJ/kg K 

r = 1.4; // Specific lieat ratio 

R = 287; // Characteristic gas constant in J/kg K 

eff_s=(rp~((r-1)/r)-1)*T01/DelT_rise; 
wf=(DelT_rise*Cp*10~3)/(u*ca*(tand(beta_l)-tand( 
beta_2))); 

disp (”%” , ef f _s*100 , ” Stage Efficiency =”); 
disp (wf,”Work done factor = ”); 


Scilab code Exa 9.10 Determination of blade and air angle 
clc ; 

u = 250; // Mean blade speed in m/ s 

rp=1.3; // Pressure ratio 

ca = 200; // Axial velocity in m/ s 

p 01 = 1 ; // Inlet pressure in bar 

T01=300; // Inlet temperature in kelvin 

Rl=0.5; // Degree of reaction 

Cp = 1.005; // Specific lieat in KJ/kg K 
r = 1.4; // Specific lieat ratio 

R = 287; // Characteristic gas constant in J/kg K 

Del_T=(rp~((r-l)/r)-1)*T01; 

//tan_betal + tan_beta2 = (R*2*u/ca); 

/ /1an_betal — tan_beta2 = (Del_T*Cp* 10" 3 / (u*ca ) ) ;. 

A=[l 1;1 -1]; B= [ (R1*2*u/ca) ;(Del_T*Cp*10~3/(u*ca)) 

] ; 

tan_beta=A\B; 

beta_l=atand (tan.beta(1)); 
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beta_2=atand (tan_beta(2)); 
alpha_1=beta_2; alpha_2=beta_l; 

disp (”degree”,beta_2,”beta2 = ”,”degree”,beta_1,” 
betal = ”); 

disp (” degree” ,alpha_2,”alpha2 = ”,”degree” ,alpha_1 , 
” alphal = ”) ; 


Scilab code Exa 9.11 Calculation of rotational speed 
clc ; 

n=4; // Number of stage 

rp=10; // Pressure ratio 

eff_p_ac=0.92; // Ploytropic efficiency of axial 

compressor 

ef f _p_cc =0.83 ; // Polytropic efficiency of 

centrifugal compressor 

Del_Trise=30; // Axial compressor stage temperature 
i n k e 1 v i n 

R = 0.5; // Reaction stage 

beta_2=20; // Outlet stator angle in degree 

D = 0.2 5 ; / / Mean d i a m e t e r o f e a c h s t a g e i n m 

wf=0.8; // Work done factor 

ca = 150; // Axial velocity in m/s 

D i = 0.3 3 ; / / I m p e 11 e r d i a m e t e r i n m 

mu=0.9; // Slip factor 

p01 = 1 . 01 ; // Ambient pressure in bar 

T01=288; // Ambient temperature in kelvin 

pif = 1.04; // Power input factor 

Cp = 1.005; // Specific heat in KJ/kg K 

r = 1.4; // Specific heat ratio 

R = 287; // Characteristic gas constant in J/kg K 

beta_l = atand (sqrt ((Cp* 10~3*Del_Trise/(wf*ca~2)) + ( 
tand(beta_2)~2))); 
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u=ca*(tand (beta_1)+tand(beta_2)); 

Nac = (u/(3.14* D)) ; 

rl=(l+n*Del_Trise/T01)~(eff_p_ac*r/(r-l)); // Total 
pressure ratio across the axial compressor 

r2=rp/rl; // Pressure ratio across centritugal 
compressor 

T02=T01*rl~((r-1)/(eff_p_ac*r)); 

T03=T02*r2~((r-1)/(eff_p_cc*r)); 

Del_Tsc=T03-T02; 

u=sqrt ((Del_Tsc*Cp*10~3)/(pif*mu)); 

Ncc=u/(3.14* Di) ; 

disp (”rps (roundoff error)”,Nac,”Speed of the 
axial compressor = ”); 

disp (”rps (roundoff er ror ) ” , Ncc , ” Speed of the 
centrifugal compressor = ”); 
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Chapter 11 

Impulse and Reaction Turbines 


Scilab code Exa 11.1 Estimation of maximum number of stages required 
clc ; 

p02=6; // Inlet pressure in bar 

T02=900; // Inlet temperature in kelvin 

p0fs = l; // Outlet pressure in bar 

eff_isenT=0.85; // insentropic efficiency of turbine 
alpha_2=75; // Nozzle outlet angle in degree 

u = 250; // Mean blade velocity in m/s 

Cp = 1 . 15* 10~3 ; // Specific lieat in J/ kg K 
r = 1.333; // Specific heat ratio 

T0fs=T02/(p02/p0fs)~((r-l)/r); // Isentropic 
temperature at the exit of the final stage 
De l_Tover all = ef f _ i senT * (T02 - T0f s ) ; // Actual overall 
temperature drop 

c2=2*u/sind (alpha_2); // absolute velocity 
c3= c2*cosd (alpha_2);// absolute vclocity 
cl = c3; // From velocity triangles 

Del_Tstage = ( c2 ~2-c 1 ~2) / (2*Cp) ; // Stage temperature 
drop 

n = Del_Toverall/Del_Tstage ; // Number of stages 
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disp (round (n),”Number of stages n =” ) ; 


Scilab code Exa 11.2 Determination of output power developed by the 
turbine shaft 

clc ; 

N=10000; // Speed of gas turbine in rpm 
T01 =700 + 273.15; // Total head temperature at nozzle 

entry in k e 1 v i n 

P01=4.5; //Total head pressirre at nozzle entry in 
bar 

P02=2.6; // Outlet pressure from nozzle in bar 

p3 = 1 . 5 ; // Pressure at trbine outlet annulus in bar 

M = 0.5; // Macli number at outlet 

alpha_2=70; // outlet nozzle angle in degrees 

D = 64; // Blade rnean diameter in cm 

m = 22.5; // Mass flow rate in kg/s 

eff_T = 0.99; // turbine mechanical efficiency 

Cp = 1.147; // Specific heat in kJ/kg K 

r = 1.33; // Specific heat ratio 

f 1 = 0.03 ; // frictional loss 

R = 284.6; // c h a r a c t e r i s t i c gas constant in J/kg K 

eff_N = l- fl; // Nozzle efficiency 
T_02=(P02/P01)~((r-l)/r)*T01; // Isentropic 
temperature after expansion 
T02 = T01 - ef f _N * (T01-T_02 ) ; // Actual temperature 
after expansion 

c2 = sqrt (2* Cp * 10 ~ 3* (T01 -T02) ) ; // Absolute velocity 

u= (3.14* D * 10 ~-2* N)/60 ; // Mean blade velocity 

// From velocity triangles 

wt2=c2*sind (alpha_2)-u; 

ca=c2*cosd (alpha_2); 

beta_2 = atand((wt2)/ca) ; 

T3=T02/(P02/p3)“((r-1)/r); // Assuming rotor losses 
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are negligible 

c3 = M*sqrt (r*R*T3); // Absolute velocity 
beta_3 = atand(u/c3) ; 
ct2=c2*sind(alpha_2); 

P=eff_T*ra*(ct2)*u/1000; // Power developed 

disp (” degree ” , beta_3 , ” Gas angle at exit = ”,”degree 
” , beta_2 , ” Gas angle at entry ” , ” ( i ) . ” ) ; 
disp (”kW (roundoff err or ) ” , P , ” Power developed = ” 
, ” ( i i ) . ” ) ; 


Scilab code Exa 11.3 Estimation of the blade angle and power produced 
clc ; 

alpha_2=65; // Nozzle discharge angle in degree 

c3 = 300; // Absolute velocity in m/s 
alpha_3=30; // in degrees 

ca2 = c3*cosd (alpha_3); // Axial velocity 
c2 = ca2 / cosd ( alpha_2 ) ; // Absolute velocity 
// ca3=ca2=ca and equal blade angles then 
ca= ca2; 

beta_2=atand((c2*sind(alpha_2)+c3*sind(alpha_3))/(2* 
ca) ) ; // Blade angle 
beta_3=beta_2; // equal blade angles 
u=c2*sind(alpha_2)-ca2*tand(beta_2); // Mean blade 

velocity 

// From velocity triangles 

ct2=c2*sind(alpha_2); 
ct3=c3*sind(alpha_3); 

WT = u* ( ct2 + ct3)/1000 ; // Work done 

sigma = u/c2; // optimum speed ratio 

eff_B=4*(sigma*sind(alpha_2)-sigma~2); 

disp (” degree” ,beta_2,” Blade angle = beta_2= beta_3 
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= ”) ; 

disp (”kJ/kg (roundoff er r or ) ” , WT , ” Power Produced 
= ”) ; 

disp (”%” , ef f _B*100 , ” Blade efficiency =”); 


Scilab code Exa 11.4 Calculation of blade angle used and the mass flow 
rate required 

clc ; 

P01=7; // Pressure at inlet in bar 

T01 =300 + 273.15; // Temperature at inlet in kelvin 

P02=3; // Pressure at outlet in bar 

alpha_2=70; // Nozzle angle in degree 

eff_N = 0.9; // Isentropic efficiency of nozzle 

WT = 75; // Power Produced in kW 

Cp = 1.15; // Specific lieat in kJ/kg K 

r = 1.33; // Specific heat ratio 

T_02=T01*(P02/P01)~((r-1)/r); // Isentropic 
temperature after expansion 
T02=T01-eff_N*(T01-T_02); // Actual temperature 
after expansion 

c2 = sqrt (2* Cp * 10 ~ 3* (T01 -T02) ) ; // Absolute velocity 
// For optimum blade speed ratio 

u=(c2*sind (alpha_2)/2); // Mean blade velocity 
beta_2=atand((c2*sind(alpha_2)-u)/(c2*cosd(alpha_2)) 

) ; // Blade angle 
// From velocity triangles 
ct2=c2*sind(alpha_2); 
w2=c2*cosd(alpha_2)/cosd(beta_2); 
w3 = w2; // Equal inlet and outlet angles 
beta_3=54; // in degrees 

ct3=w3*sind(beta_3)-u; 

m=(WT*10"3)/(u*(ct2 + ct3)); // Gas mass flow rate 
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disp (” degree ” ,beta_2,” Blade angle = ”); 
disp (” kg/s” ,m, ” Gas Mass Flow Rate = ”); 


Scilab code Exa 11.5 Determination of gas temperature velocities and dis- 
charge angles at the blade root and tip radii 

clc ; 

P01=4.6; // Total head inlet pressure in bar 

T01 =700 + 273.15; // Total liead inlet temperature in 

k e 1 v i n 

P2=1.6; // Static head pressure at mean radius in 

bar 

Dra_h = 10; // Mean blade diameter / blade heiglit 

1 c = 0.1; // Nozzle losses coefficient 

alpha_2=60; // Nozzle outlet angle in degree 

Cp = 1.147; // Specific lieat in kJ/kg K 
r = 1.33; // Specific lieat ratio 
ra = 20; // Mass flow rate in kg/s 

R = 284.6; // c h a r a c t e r i s t i c gas constant in J/kg K 

T_2=T01*(P2/P01)~((r-1)/r); // Isentropic 
temperature after expansion 
T2 = (lc*T01+T_2)/(1 + lc) ; // Actual temperature after 
e x p a n s i o n 

c2 = sqrt (2*Cp*10~3* (T01-T2) ) ; // Absolute velocity 
// From velocity triangles 
ca=c2*cosd(alpha_2); 

row = P2 * 10~5/(R*T2) ; // Density of gas 
A=m/(ca*row); // Area 

Dm=sqrt (A*Dm_h/3.14) ; // Mean Diameter 

h = Dm/10; // Blade heiglit 

rm = Dm/2; // Mean radius 

// At root 

r_root=(Dm-h)/2; 

/ / At t h e t i p 
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r_tip = (Dra + h)/2; 

// Free vorte flow 

ct_mean=c2*sind (alpha_2); 

// At the root 

ct2_root = (ct_mean*rm)/r_root ; 
alpha2_root=atand(ct2_root/ca); 
c2_root=ct2_root/sind (alpha2_root); 
T2_root=T01-c2_root~2/(2*Cp*10~3) ; 

// At the tip 
ct2_tip=ct_mean*rm/r_tip; 
alpha2_tip = atand (ct2_tip/ca); 
c2_tip=ct2_tip/sind(alpha2_tip); 

T2_tip = T01-c2_tip“2/(2*Cp*10~3) ; 

disp (”degree” ,alpha2_root ,” Discharge angle at the 
root = ” , ”m/s” , c2_root , ” Gas velocity at the root 
= ” , ”K” , T2_root , ” Gas Temperature at the root = ” 
”A the Root”); 

disp (” degree”,alpha2_tip, ” Discharge angle at the 
tip = ” , ”m/s” , c2_tip , ” Gas velocity at the tip = 

, ”K” , T2_tip , ” Gas Temperature at the tip = ”,”A 
t h e t i p ” ) ; 
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